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Chapitre 1 

Introduction: les matrices aleatoires en 
physique et en mathematiques 

Je vais presenter dans cette habilitation, des resultats obtenus ces dernieres annees sur le modele 
a deux matrices et les families de polynomes biorthogonales. 

0.1 Bref historique du sujet 

Les modeles de matrices aleatoires ont ete introduits en 1951 par E. Wigner [7H] . en physique 
nucleaire, comme modele pour la repartition des niveaux d'energie des noyeaux lourds. Cette 
theorie a eu un succes considerable, et s'est imposee comme outil tres puissant pour etudier 
de nombreux autres phenomenes physiques [H3 |2TJ HUB HI] • En vrac: chaos quantique, conduc- 
teurs mesoscopiques, croissance des cristaux, problemes de frontieres libres (probleme de Saffman- 
Taylor), gravitation quantique, theorie des cordes, repliement des proteines, etc... La raison de ce 
succes tient dans les proprietes d'integrabilite des modeles de matrices. Les matrices aleatoires 
ont joue un role important aussi en mathematiques [SHI El ESI- Elles sont etroitement reliees aux 
polynomes orthogonaux, aux problemes isomonodromiques, probleme de Riemann-Hilbert. Les 
integrates de matrices sont le prototype des fonctions r. De plus, recement, on a compris que 
toutes ces notions avaient egalement beaucoup a voir avec la geometrie algebrique. 

Le modele a une matrice (2H EH] , qui est relie aux polynomes orthogonaux habituels, a ete 
beaucoup etudie. Toutefois, il ne donne acces qu'aux courbes hyperelliptiques, et aux systemes 
d'equations different ielles d'ordre 2. II aparaissait interessant et important de generaliser les 
notions qui aparaissaient dans le modele a une matrice, a un cas un peu plus general. C'est 
pourquoi le modele a deux matrices a alors attire l'attention des physiciens et des mathematiciens. 

Le modele a deux matrices, est relie a deux families de polynomes biorthogonales. Ces 
deux families satisfont deux systemes differentiels d'ordre differents, dont on peut montrer qu'ils 
sont duaux l'un de l'autre. La courbe spectrale de ces systemes, est une courbe algebrique, pas 
necessairement hyperelliptique, en particulier, elle peut avoir des singularites rationelles y ~ x p l q 
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avec (p, q) entiers quelconques. 

Du point de vue des physiciens, le modele a deux matrices a ete introduit pour modeliser le 
modele d'Ising [221 sur un reseau aleatoire, i.e. en termes de theories de champs conformes, un 
modele minimal (3,4). II a rapidement ete compris ^H] que le modele a deux matrices donnait 
acces a tous les modeles minimaux rationels (p, q). C'est dans ce cadre qu'ont ete inventees et 
etudiees les "equations de boucles" i(T0j [2Tj . et le developpement topologique |B5j . Plus recement, 
les modeles de matrices ont joue un role tres important en theorie des cordes |2"2"] . 

0.2 Principaux resultats obtenus 

Parmis les resultats les plus importants que j'ai obtenu, avec mes collaborateurs, sur ce sujet, et 
qui sont presentes dans cette these: 

• L'expression des fonctions de correlations de valeurs propres en terme de determinants de noyaux 
de polynomes bi-orthogonaux ( theoreme HUH. 11 [P10]). 

• Un theoreme de Christoffel-Darboux pour ces noyaux ( theoreme 11114.21 [P8]). 

• L'obtention de deux systemes differentiels pour les deux families bi-orthogonales de 
polynomes, leur expression explicite et le calcul de leur trace ( theoreme III IS . 11 theoreme IHJ6.11 
theoreme IIII6.21 [P2]), et un premier pas vers le calcul de la fonction r isomonodromique [P7]. 

• La dualite spectrale entre les deux systemes differentiels ( theoreme IIII7.21 [P8]). 

• L'expression exacte de la solution fondamentale de ces systemes differentiels ( theoreme IIH8.11 
[P6]). 

• Les asymptotiques a x grand des solutions fondamentales ( theoreme IIH9.11 [P6]), et les matrices 
de Stokes associees. 

• La propriete de dualite spectrale a permis d'ecrire un probleme de Riemann-Hilbert associe 
(paragraphe IIlHUJ [P6]). 

• L'expression d'une fonction de correlation mixte (i.e. qui ne peut pas s'ecrire directement en 
termes de valeurs propres des deux matrices) ( theoreme III112.11 [P4]). 

• La limite iV grand des integrales de matrices formelles, en particulier l'obtention d'une equation 
algebrique pour la resolvante ( theoreme IVI1.31 theoreme IV12.11 [Pll], [P9], [P3]), et le calcul 
de certaines fonctions de correlations mixtes dans la limite N grand ( theoreme IVU.41 [Pll], 
[P9], [P3]). 

• Le developpement a N grand des integrales de matrices formelles ([Pll], [P9], [P3], [P5], 
[PI]) en particulier, le fait que le terme enl/iV 2 du developpement de l'energie libre est relie au 
determinant d'un Laplacien sur une courbe algebrique ( theoreme IVI3.11 theoreme IVI3.21 [P5], 

[PI])- 

• Une conjecture, (tout de meme basee sur de serieux arguments heuristiques) pour les asympto- 
tiques a N grands des polynomes biorthogonaux (conjecture V 12.ll [Pll], [P9]), du type Deift & 

co [U]. 
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0.3 Plan 



La these est organisee comme suit: 

- Le chapitre 2 donne les definitions du modele a deux matrices, et explique les quantites que 
Ton cherche a calculer. 

- le chapitre 3 presente la methode des polynomes biorthogonaux, et les systemes differentiels 
associes. Un certains nombres de resultats comme la dualite, les asymptotiques, y sont exposes. 

- le chapitre 4 presente la methode des equations de boucles, et la courbe algebrique associee. 
Un certains nombres de mes contributions comme le calcul du developpement en 1/N 2 de l'energie 
libre y sont exposes. 

- le chapitre 5 montre, comment en combinant les deux methodes precedentes, on peut for- 
muler une conjecture pour les asymptotiques iV-grand des polynomes biorthogonaux, en termes 
de geometrie algebrique. 

- le chapitre 6 est une conclusion qui propose les suites possibles envisageables de ces travaux. 
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Chapitre 2 



Trois definitions du modele a deux 
matrices 



II existe plusieurs definitions du modele a deux matrices, non equivalentes. Suivant le contexte, 
c'est a dire a quelle application des matrices aleatoires on s'interesse, il faut choisir l'une des 
definitions. Dans la limite de grande taille des matrices, a l'ordre dominant, les trois definitions 
conduisent souvent aux memes resultats, mais ce n'est pas vrai au dela de l'ordre dominant. De 
plus un certains nombre de methodes de calcul sont communes aux trois modeles, et les trois 
modeles font appels aux memes concepts de geometrie et integrabilite. De la, il sort que les trois 
modeles ne sont pas toujours bien distingues dans la litterature, et Ton rencontre de nombreuses 
confusions. 

Nous allons introduire ces trois modeles, et expliquer comment ils sont relies entre eux (le 
modele 1 et le modele 2 sont relies par une transformed de Fourrier, et le modele 2 et le modele 3, 
par leur developpement en serie asymptotique). 

Avant de les introduire, il faut rappeler quelques notions de base. 

1 Quelques definitions utiles 

1.1 Mesure de Lebesgue pour les matrices hermit iennes 

On definit la mesure sur Hn (=espace des matrices hermitiennes de taille N): 



M e H 



N 



dM := — 





ou Un est le "volume du groupe unitaire": 




(II.1-2) 
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1.2 Changement de variable: matrice valeurs propres + unitaire 

Toute matrice hermitienne M G Hn peut s'ecrire J40J (pas de fagon unique, i.e. permutation des 
valeurs propres Sjv et phases diagonales U{1) N ): 

M = C/AC/ t (II.1-3) 

ou A = diag(Ai, . . . , A at) est diagonale, et U G t/"(-/V) est unitaire. On a: 

dM = A(A) 2 d£/dA (II. 1-4) 

oil dA = rijdAj, dU est la mesure de Haar normalisee sur U(N), A(A) est le determinant de 
Vandermonde: 

A(A) := - A,) (II.1-5) 

i>j 

Pour A/" = 1, on definit A := 1. 

1.3 Integrates sur le groupe unitaire: formules de Itzykson-Zuber- 
Harish Chandra, Morozov et Eynard Prats Ferrer 

Soient A = diag(ai, . . . , on) et B — diag(&i, . . . , 6jv) deux matrices diagonales non degenerees, et 
dU la mesure de Haar sur le groupe U(N), on a: 

Theoreme 1.1 Theoreme (Eynard, Prats-Ferrer [36]): 

Pour toute fonction polynomiale invariante ( i. e. produit de traces de produits des deux matri- 
ces) F, on a: 



[ F(A,UBU^)e- T * AUBin dU 

Ju(N) 



'U(N) 

J2 (_i)- e -Tr^ f e -Tvrrt F (A T + TjBa + Tl)dT (II.1-6) 



NlA(a)A(b) c 

oil T(N) est V ensemble des matrices complexes triangulaires superieures strides, muni de la 
mesure produit des mesures de Lebesgues des parties reelles et imaginaires de tous les elements de 
matrice. Si a est une permutation, B a est la matrice diag(6 CT (j)), et la constante vaut: 

cn = - (in-?) 

En particulier avec F = 1 on a: 
Theoreme 1.2 Formule Itzykson-Zuber -Harish Chandra \4 1\ \4 b\ \T®i : 

u(N) A(o)A(6) 
oil la matrice E est donnee par: 

E^ :=e a ^ . (II. 1-9) 
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Et, avec F(A, B) = Tr — ^ -^-^ vue comme serie formelle en puissances de 1/x et 1/y: 

Theoreme 1.3 Theoreme de Morozov simplifie (eynard [29]): 

La fonction generatrice definie par le developpement formel en puissances de 1/x et 1/y: 

tr _ L^_Lm\ := y _J 1 /\ a . |2\ (ii. no) 



est donnee par: 



(\u--\ 2 \ ■= Ju{N) — aii-m 

JU(N) UU e 



«jhi u vh W ) = 1 - det ( 1 ^^^ B " 1 (IU " 12) 



Les correlations (\Uij\ 2 ) ont ete calculees par A. Morozov jM], puis simplifiee par moi meme dans 
[2H], et utilisee pour le resultat obtenu par Bertola-Eynard dans [P4] et presente au paragraphe 
H2]dans le cadre de cette habilitation. 

1.4 Ensemble des matrices normales sur un chemin 

Definition 1.1 Pour tout contour 7 du plan complexe, on definit Vensemble de matrices H n {^f): 
# n ( 7 ) := {M G GL n (C) 1 3A = diag(Ai, . . . , A n ) , A» G 7 , 3Z7 G U(n) , M = C/A^} (II.1-13) 
On le munit d'une mesure analogue a eg. \II.l-4 ): 



dM = A(\ydUdA (II.1-14) 

L'ensemble des matrices hermitiennes est H„ = H n (R). 
Notons que si M G H n (j), on a: 

[M,M f ] = (II.1-15) 

d'ou le nom de matrices normales [7^j . 

2 Le modele hermitien et le modele normal 

Le modele hermitien est relie a deux families de polynomes biorthogonaux (que Ton etudiera 
dans le chapitre EJ) ■ Une generalisation naturelle des polynomes biorthogonaux (dits polynomes 
semi-classiques, cf [3]), nous conduit a introduire une generalisation du modele hermitien, que 
j'appelerai le modele normal sur un contour. 
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2.1 Les donnees 



Donnons nous deux polynomes V% et V2 de degres respectifs di+1 et c?2+l, a coefficients complexes, 
appeles potentiels: 

di+l d2+l ~ 

^(x) ^o + Ef^. ^(y) = a>+£fy fc , ( IL2 - 1 ) 

fc=i fe=i 

un entier iV > min (c?i, c/2), et un nombre complexe non nul T, appele temperature. 

II existe d\ chemins homologiquement independants allant de 00 a 00, T x k, k — 1, . . . , d\, tels 
que l'integrale: 

e~^ Vl ^ dx (II.2-2) 

soit absolument convergente. II suffit que chaque chemin parte de 00 dans un secteur ou 
0, et revienne a 00 dans un autre secteur ou Re^^ > 0. 

De meme, il existe d 2 chemins homologiquement independants allant de 00 a 00, T y k, k = 
1, . . . , c?2) tels que l'integrale: 

f e -^v 2 (y) dy (II.2-3) 
Jvy k 

converge. II suffit que chaque chemin parte de 00 dans un secteur ou Re^ > 0, et revienne a 
00 dans un autre secteur oil Re^r^- > 0. 

Deux telles bases de chemins etant choisies, donnons nous une matrice Kij a coefficients com- 
plexes, non identiquement nulle, et definissons: 

r := X>j f i x (II.2-4) 

hi 

Remarque 2.1 nous aurions pu considerer un modele plus general, ou V\ et V2 ne sont pas des 
polynomes, mais tels que V[ et V{ soient des fractions rationelles. Dans ce cas, il faudrait introduire 
des bases d'homologies de tous les contours sur lesquels J dxe~~ Vl ^ a un sens. Cela inclut les chemins 
allant d'un pole de V\ a un autre pole, et cela inclut aussi des segments finis,... C'est le cas semi-classique 
introduit par Nous ne le discuterons pas ici par souci de clarte et simplicite. Le cas general est 
presente dans jZj, il n'est pas tres different du cas polynomial presente ici. 

2.2 Modele hermitien 

C'est le cas oil il est possible de choisir les bases d'homologies telles que: 

T = R x R (II.2-5) 

c'est a dire si et seulement si les conditions suivantes sont satisfaites: 

d\ et d 2 sont impairs, et g^+i/T et gd 2 +i/T ont une partie reelle strictement positive. Et dans 

le cas ou di = d% = 1, la forme quadratique h { ^ ~ ) a des valeurs propres de parties reelles 

V 1 92 J 
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strictement positives. Autrement dit si et seulement si Re [^(Vi(x) + V 2 (y) — xy)^j est bornee 
infer ieurement sur R x R. 

On introduit alors 

Definition 2.1 Mesure du modele hermitien sur x H^: 

d f i(M 1 ,M 2 ) := ^J— e -^mMi)+V2(M a )- Ml M a ) dMidM2 (IL2 _ g) 

■^Herm 

oil la normalisation Zn crm 



^Herm == e"f^ H - := / e"T tr mM 1 )+V 2 (M,)-M,M 2 ] ^ ^ ( jj 2 _ ?) 

J H N xH N 

est appelee fonction de partition, et F^ eim est appelee energie libre. 
En ecrivant comme eq. (JII.l-3j) : 

M 1 = UXU ] , M 2 = [UV) Y (UVy (II.2-8) 
avec X = diag(xi, . . . , xn) et Y = diag(yi, . . . , Un)i et U et V unitaires, on a: 

d»(M u M 2 ) = lA{x) 2 A(y) 2 e-^ tr ^ +v ^e^ XVYV 'dXdYdUdV (II.2-9) 
Z 

En utilisant |H.l-8l l'integrale sur U et V donne une mesure induite pour les valeurs propres Xi de 
Mi et yi de M 2 : 

Definition 2.2 Mesure sur K N x R^ pour les valeurs propres du modele hermitien: 

N N 

'■±{x)^{y) -| 

N\ ^Herm 



i=l i=l 

La normalisation 

N 



Znerm - 4 / det (eT^i) A(2r) A(j/) TT e"^ W*')+ v Md Xi d Vi 

N 

A(x) A(y) We-^^^^-^dxidyi 

:= e"^ FHcrm (II.2-11) 
est appelee fonction de partition, et Fnerm est appelee energie libre. 
Notons que 



N v N(N-l)/1 

Zuerm = -^Herm ( ) (II. 2-12) 
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Remarque 2.2 L'espace des parametres du modele hermitien est de dimension d\ + cfo + 4: les d\ + 2 
coefficients de Vi, les ofo + 2 coefficients de V2, la temperature T, et il faut retrancher la redondance des 
coefficients constant de V\ et V2. 

Remarque 2.3 Le modele de matrices hermitiennes, est celui qui peut s'appliquer a des problemes de 
physique de la matiere condensee j33 EE1 E01 • En effet les matrices qui apparaissent dans ces problemes 
representent des Hamiltoniens, et sont toujours hermitiennes (elles peuvent avoir aussi des symetries 
suplementaires) . 

2.3 Le modele Normal sur une classe d'homologie de contours T 

II apparait naturel d'etendre la definition de la fonction de partition Zn eTm , donnee par la seconde 
partie de l'equation eq. (|II. 2-11)1 a n'importe quel contour T comme defini par eq. f)II.2-4j) . 
On definit done, comme generalisation de eq. (jll.2-10|) . la mesure pour les valeurs propres: 

Definition 2.3 Mesure sur YliLi ^fo x T y i t pour les valeurs propres du modele normal sur un 
chemin T: 



C^Norm, ki,...,k N ;h,...,l N 

(x u ...,x N ;y!,.. .,y N ) 
A(x) A(y) det ( K ki ,ij e ^ Xi 



N\Z Norm (n) 



I J e -- [Vl{x * )+VM dx t dy t (II.2-13) 



ou ZNorm(^) = e t 7 Norml ' est la fonction de partition, et F^ orrn {K) est Venergie libre. 
On a: 

:= J_ e -f^N 0rm M 

N „ N 



ki,...,k N ,h,-,lN 1=1 JXl ^ 1 k i j=l ' 

A{x) A{y) det (/^.e^*) f[ e^ 1 ^^ 1 d Xi d Vl 

i 

A(x) A(y) Y[K kiM e^i J]e"* [*(*<)+v.(i«)] dXidyi 



{II. 2 - 14) 



ou Ton a decompose le determinant comme une somme sur les permutations de SV- L'integrale 
a effectuer est la meme pour chaque permutation (il suffit de renommer y ai — > yi et l ai — > k), et 
done: 

AT 



Z N orm(«0 = II II / 



fci fcwJi z w i i=i Jx i^k^vy h 
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(//.2 - 15) 

Remarque 2.4 Notons que d'autres generalisations du modele hermitien sont possibles, en particulier, 
on pourrait choisir de mettre des valeurs absolues aux Vandermonde carres. La definition ci-dessus est 
celle qui s'interprete en termes de polynomes biorthogonaux presentes au chapitre C3 et surlaquelle ont 
porte une grande partie de mes travaux. 

Remarque 2.5 En rassemblant les paires (xi,yi) qui sont sur le meme produit de chemins (T x k x T y i), 
on peut ecrire: 

w«) - e n^n<rnn/ / 

£ M n k , l= N 1 lk > 1 k,i k,i i=i ■ / *fc,i,ier** Vwer*i 

A(s) A(y) J] JJ e -^W(«w.O+^(v fc .w)-*w.iW*,w]d a . Jfc)W d y m 

k,l 4=1 

rife ; 

e Tr^n<rnn/ / 

Efeji „ m =jv 1 lk > 1 k ' 1 ' k,i k,i i=i • / ^.-' ' '• ■ / w*,«,ier» I 

A(x) A(y) H [] e-fW^H^w)]^ d y h>l>i 

nk > 1 ' k,l 4=1 

(JJ.2 - 16) 
et Ton reconnait des integrales sur C/ nfc ; : 



n k ,i 



= n\ n — n^^.-D / / 

e -f tr [^i(M., fc> ,)+V a (M 1/>fc , I )-Af 1 ,, fc , I W I ,, fc , ! ] dMjBjfc)J 

]~[ det (M x>k ,i ® ln fc , ; , - ln fe ,( ® M Xjk > ; i>) 
(k,l)>(k',V) 

Yl det (Mj, )fe)Z (8) 1^, - <g> M y>k r )V ) 

(k,l)>(k',V) 

{II. 2 - 17) 

ou Ton a ordonne les paires (k, I) par ordre lexicographique. 



Remarque: L'espace des modules du modele normal est de dimension d\ + d 2 + 3 + d\d 2 . Les 
modules sont les d\ + 2 coefficients de Vi, les d 2 + 2 coefficients de V 2 , la temperature T, et les 
did 2 coefficients K k> i, et il faut retrancher la redondance des coefficients constant de V\ et V 2 et un 
facteur global pour k. 
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3 Le modele Normal a symetrie brisee 



Nous allons considerer individuellement chaque terme de la somme eq. (jll.2-17|) . autrement dit on 
se place a nombre de valeurs propres fixe sur chaque contour. Cela revient a briser le groupe de 
symmetrie U(N) en Y[ k i U(n k ,i)- Comme nous le verrons au chapitreEl ce modele est tres utile 
pour faire le lien entre polynomes biorthogonaux et developpement a N grand. 

3.1 Les donnees 

On se donne deux entiers d\ et d 2 superieurs ou egaux a 1, et un entier N > min (di, d 2 ). 
On se donne d\ x d 2 entiers n k[ , 1 < k < d±, 1 < / < d 2 tels que 



J2n k)i = N (II.3-1) 



k,l 

et, pour des raisons qui deviendront claires au chapitreEl on note: 

e k,i := ^ (H.3-2) 

Les e k j sont appeles fractions de remplissage. On note 

n := {n k) i} et e := {e k ,i} (II. 3-3) 

Donnons nous aussi deux polynomes V\ et V 2 de degres respectifs d\ + 1 et d 2 + 1, a coefficients 
complexes, appeles potentiels: 

di+l d 2 +l ~ 

v x (x) = 9 o+J2j xk > V M = 9° + E f y k ' ( IL3 - 4 ) 

k=l k=l 

et un nombre complexe non nul T appele temperature. 

Et on se donne des bases de classes d'homologies de chemins dans C definies comme au para- 
graphe precedent: 

f fe k = l,...,di , f/Z = l,...,d2 (II.3-5) 

3.2 Le modele 

Definition 3.1 Mesure du modele normal d symetrie brisee sur Ylk 1 Hn k XT* k) x H nkl (T y i): 

Y[ det (M xA i ® l nfc , ; , - l nw <g> M XtWJi ,) 

(k,l)<(k' ,/') 

J] det (M„ ffe) , ® 1^, - l nfcl , ® ^i/,fc',/0 
(jfe,0<(fc',J') 
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(//.3 - 6) 

oil les paires (k, I) sont ordonnees par ordre lexicographique, et ou la fonction de partition et 
I'energie libre sont donnees par: 



^Normb(^) '■— £ T 



V ^Normb(«) 



(//.3 - 7) 



-Q e -f tr mM*M)+V>(M y M)-M a MM v , k , l ] dMxM dMy k i 
k,l 

Yl det (M X)k j ® l ny i , - l nk l <g> M X)k r t ii) 

(k,l)<{k',l') 

Yl det (M yAl ® l nfc , it< - l nfeji ® M Vifc /,i') 

(k,l)<{k',l>) 



En utilisant la decomposition eq. (jll.l-3|) . et en integrant sur les groupes unitaires avec la formule 
eq. pi. 1-8(1 . on obtient la mesure induite pour les valeurs propres 

Definition 3.2 Mesure pour les valeurs propres du modele normal a symetrie brisee, sur 
Uk,i (r x fc x T y i) nk ' 1 (i.e. Vfc = l,...,d x , VZ = 1, . . .,d 2 , Vi = 1, . . . ,n k ,i , x k ^i e T x k , Vk,i,i e 



ZMnr m h W - LJ - \ V / AA 



'NormbW , =1 



(//.3 - 8) 

cm /a fonction de partition et I'energie libre sont donnees par: 

f n ^ ■— g— ^r-PNormbC"') 

n k,l 

■= II II 



J NormM 



A(ar) A(y) JJ [J e~T I^m^+^m.O-^,^,^^. ^ 



Jk,U 

k,l i=l 

(J/.3-9) 

3.3 Relation avec le modele Normal 

II est clair que le modele normal est tel que: 



Z N orm(«0 = N\ ( II K Z" J ^NormbM (II.3-10) 

n \ k,l J 
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Autrement dit, les k^i sont les variables duales des nk,u e t les deux modeles sont relies par une 
transformed de Fourrier (avec la contrainte ^2 n k,i — N) . 

Remarque: L'espace des modules du modele normal brise est de dimension d\ + d 2 + 3 + d\d 2 , 
comme le modele normal. Les d\d 2 — 1 coefficients Kk,i (moins la normalisation redondante) ont 
etes remp laces par les d\d 2 — 1 fractions de remplissages independantes e^. D'apres eq. (jll.3-10|) . 
les Kk,i apparaissent comme des variables duales des e^i- 

4 Le modele formel 

Le modele formel pourrait aussi etre appele modele combinatoire. En effet, il n'est pas defini par 
une integrate de matrices (ni de valeurs propres), mais par une serie formelle en puissances de la 
temperature T, qui sert de serie generatrice pour le denombrement de certains types de graphes 
[2TJ EH UZ1 123 SJ] ■ Pour ce modele, les questions de convergences et de bases d'homologies ne se 
posent pas. 



4.1 Les donnees 



On se donne deux entiers d\ et d 2 superieurs ou egaux a 1. Donnons nous aussi deux polynomes 
V\ et V 2 de degres respectifs d\ + 1 et d 2 + 1, a coefficients complexes appeles potentiels: 



V 1 (x)=go+Y2jX k 

k=l 



V 2 (y) 



d 2 +l ~ 

9k k 



k=l 



(II.4-1) 



et deux nombres complexes non nuls T et N. T est appele temperature. 

Le potentiel V\(x) + V 2 (y) — xy possede did 2 extrema dans C x C, que Ton note: (xj,yi 
I = 1, . . . , d\d 2) tels que: 

' V{(xj) = Vl 
V 2 \yi) = x I 

et on se donne d\d 2 nombres complexes arbitraires ej , I = 1, . 
de remplissages, tels que: 



(II.4-2) 

d\d 2l que l'on appelle fractions 

(II.4-3) 



i=i 



Definition 4.1 on definit les constantes de couplages: 



VI = 1,.. 


. , d\d 2 


yk = 


2,.. 


,di + l 


9k,i 




VI = 1,.. 


. , d\d 2 


\/k = 


2,.. 


,d 2 + l 


9k,i 




Vz' j 




yk = 


1,.. 


,oo, 


hk,i 




Vi ^ j 




yk = 


1,.. 


,oo, 


hk,i 




Vz j 




yk,i 


= 1,- 


..,oo, 


hk,I;l,J 




Vz' ^ j 




yk,i 


= 1,- 


..,oo, 


hk,I;l,J 





fc-1! 
fe-1! 

Li Hi (yj-yi) k 

fe+Z-l! 1 



(II.4-4) 



fe— 111— l! (xj—xj) k (xj—xj) 1 

fc+l-1! 1 

fe— 1! /— l! (yj-yi) k (yi-yj) 1 
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Notons que: 




(II.4-5) 



Remarque 4.1 Notons, que contrairement au modele precedent, N et nj = Nei ne sont pas 
necessairement des entiers. 

4.2 Combinatoire de surfaces discretisees 

Definition 4.2 SoitQ V 'ensemble des graphes fermes (sans bords), pas necessairement connexes, 
formes de polygones orientes (une face superieure et une face inferieure), colles le long de leur 
bord ou par leurs centres, avec les regies suivantes: Chaque polygone a k cotes (k > 1 on peut avoir 
des 1-gones et des 2-gones) porte un "signe" (+ ou —) et une "couleur" (1 = 1,... ,did 2 ). Deux 
polygones peuvent etre colles par leur bord seulement si Us ont la meme couleur, ( et en respectant 
V orientation des faces). Deux polygones peuvent etre colles par leurs centres seulement si Us sont 
de couleur differente et de meme signe. 

Definition 4.3 Pour un graphe G G Q, on note: 

• nkj(G) = nombre de k—gones de signe + et de couleur I. 

• nk,i(G) = nombre de k—gones de signe — et de couleur I. 

• n p j(G) = nombre d'arretes communes a deux polygones de couleur I et de signes quelconques. 

• n ++ j(G) = nombre de paire de polygones de signe + et +, de couleur I colles par un bord. 

• n ,i(G) = nombre de paire de polygones de signe — et — , de couleur I colles par un bord. 

• n<k,i;i,j(G) = nombre de paire de polygones formees d'un k—gone de signe + et de couleur I , et 
d'un I— gone de signe + et de couleur J, colles par leurs centres. 

• nk,r,i,j(G) = nombre de paire de polygones formees d'un k—gone de signe — et de couleur I , et 
d'un I— gone de signe — et de couleur J, colles par leurs centres. 

• li(G) = nombre de sommets de couleur I. 

• #Aut(G) = cardinal du groupe des automorphismes de G. 

• x(G) = Caracteristique d'Euler-Poincare de G, en comptant les liens entre centres de polygones 
comme des cylindres. x(G) est toujours pair. 




(II.4-6) 



k 



i 



• et nr{G) est donne par: 




(II.4-7) 



KJ k>l l>l 
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Remarque 4.2 On a l'inegalite: 

n T (G) > i (X)E( n M + »M) + Z)Z)Z)( n V;^ + ft *,-r;i,j)) ( IL4 " 8 ) 

\ 7" fc=l KJ k>l 1>1 J 

Definition 4.4 Le pozrfs Feynmann d'un graphe G est: 

:= #AutfG) T " T(G) n n (»w + ?x,r-' (G) n + Th kJ r AG) 

^ ' 7=1 fc=3 fc=3 

di<72 <ii(fe 

n c/ (g x 2 / (g) n c/ (G) n^/ (G) ^?/ (G) n ^*/ (G) 

7=1 fc>di+l 7=1 fc>d 2 +l 

oo oo 

n n n «y (G) ^y (G) n e " (G) n ^ ++/ ^ - i r npAG) 

I<Jk=ll=l I I 

(HA - 9) 
4.3 Le modele 

On definit la fonction de partition du modele formel comme une serie generatrice (au sens combi- 
natoire), i.e. une serie formelle en puissances de T: 

Definition 4.5 Fonction de partition du modele formel 

Z Faim (e) := Y\e-^ tl(Vl{xi)+V2(yi) ~ XlVl) Hixi-Xj^'Ylfa-yj)"*™ 

I KJ KJ 

^W(G) (II.4-10) 

Gag 

Notons que l'inegalite III.4-8I implique que pour chaque puissance de T, seul un nombre fini de 
graphes contribuent, et done III. 4-101 definit bien une serie formelle en puissances de T: 



oo 



Z Form = Y,A n T n (II.4-11) 



n=0 



On definit aussi l'energie libre formelle: 
T 2 

-^Form • = jy2 ^ ■^F° rm 



^2 Te i i v ii x i) + v 2(yi) - xiyi) -t 2 ^2 e i e J ln ( x i - x j) - t 2 z2 e/6j ln & ~ yj "> 

I 

n 

N 2 



I KJ KJ 

w £ w ^ ( IL4 - 12 ) 

conn 



qui est aussi une serie formelle en puissances de T, qui est la meme somme que pour Z Form , mais 
reduite aux diagrammes connexes seulement: 



F FoTm = J2B n T n (II.4-13) 



n=0 
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4.4 Relation avec le modele normal a symetrie brisee 

Ces definitions formelles viennent de l'interpretation combinatoire de II 1 .3-71 avec nk = Nek, par la 
methode des diagrammes de Feynman, comme serie generatrice de diagrammes (ces graphes sont 
introduits dans [13] )■ En effet, on reecrit eq. ()11.3-7|) sous la forme: 



^Normb 



/d\d,2 
1=1 

Yl det (M x j ® 1 ®M x ,j) det(M y>/ ®l nj In j 



{II A - 14) 



KJ 

d\d2 

\\ dM xJ dM y> i e~ s 
i=i 



puis on pose M x j = xjl ni + X T , M y j = yj\ ni + Yj, et Taction S, ecrite en termes des matrices 
X[ et Yi est: 



^(«|itr.Yf + i|itrlf-trAw) 



1 \fc>l fc>l 



7"<J fe>l i>l 



A-/ 



+£££% £ «> 1 'J 



(J/.4 - 15) 



7<J k>l 1>1 



La methode de Feynmann consiste a developper en serie de Taylor l'exponentielle de tous les 
termes sauf la partie quadratique de la premiere ligne, puis a calculer les integrales gaussiennes 
restantes par le theoreme de Wick (Ml 13 El EU EH| ■ Les chemins d'integration ne jouent aucun 
role ici. Seuls les voisinages des points xi, yi jouent un role. La methode de Feynman associe 
a chaque terme du developpement, un diagramme G G Q, avec un poids qui est celui ecrit plus 
haut. Ceci est presente en appendice de [T5] . 

Toutefois, il faut savoir que, meme lorsque N et les ni = Nej sont des entiers, les definitions 
111.4-101 et 111.3-71 ne coincident pas toujours. 

Elles ne coincident que dans certains cas (pour des ensembles ouverts de potentiels V±, V2, pour 
certaines valeurs de fractions de remplissage, et pour des choix appropries de contours T x k et T y i)- 
II faut en particulier que Ton puisse trouver une base de chemins T x k,T y i telle que chaque paire 



19 



de chemins passe par un seul extremum (xi, yi) et telle que le potentiel effectif pour une paire de 
valeur propre de (Mi,M 2 ) ait un unique minimum. 

Lorsqu'elles coincident. III.4-101 est le developpement limite de III3-7I en puissances de T. En 
general, ce developpement est un developpement asymptotique, qui ne converge pas vers une 
fonction unique [2TJ I2Z] • 

4.5 Developpement topologique 

Pour chaque puissance de T, la serie 111.4-131 compte un nombre fini de diagrammes connexes, et 
done chaque B n est un polynome en 1/N 2 : 

deg B„ 

B n (N) = J2 B n , h N- 2h (II.4-16) 

fc=0 

Ceci autorise a definir les series a puissances de N fixee (i.e. a topologie fixee puisque la puissance 
de N est la caracteristique d'Euler). 

oo 

F£L--=J2 B ^ Tn ( IL4 - 17 ) 

n=0 

Chaque F^' est une serie formelle en puissances de T, qui compte les diagrammes de Q avec le 
meme poids que pour Z Form , restreints aux diagrammes connexes de genre h OH EDI- Ainsi, 
par definition du modele, on a le developpement de l'energie libre en puissances de N~ 2 : 

Definition 4.6 Developpement topologique de l'energie libre du modele formel: 



F Form = J2 F £LN- 2h (II.4-18) 



L p(h) AT -2h 
h=0 

II s'avere que chaque F^l m est une fonction algebrique de tous ses parametres (les coefficients de 
V\ et V 2 , les ek, T), et done possede un rayon de covergence non nul comme serie en puissances 
de T. Une partie de mon travail presente pour cette habilitation, a consiste a calculer -Fporm' e ^ ^ 
l'identifier avec le determinant du Laplacien sur une courbe algebrique [P5], [PI]. Au chapitrelU 
nous verrons l'expression de Fporm e ^ -^Form- 

Par contre, la somme eq. (jll.4-18|) est une serie asymptotique en puissances de 1/N 2 , et, en 
general, ne se resomme pas en une fonction analytique des parametres. 

Remarque: L'espace des modules du modele formel est de dimension d\ + d 2 + 3 + did 2 , comme 
le modele normal et le modele normal brise. 

Remarque: Le modele formel est celui qui est utilise pour les applications des matrices aleatoires 
au denombrement des surfaces discretisees, et subsequement, pour les applications aux theories 
des champs conformes, a la gravitation quantique et a la theorie des cordes en physique [2H EH 

nzumEniEa. 
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5 Observables 



Nous allons maintenant decrire les quantites que l'on cherche a calculer dans ces modeles. Lorsque 
cela est possible, on cherche bien sur a calculer ces quantites pour tout N, mais l'objectif principal 
en physique est de les calculer pour N grand, dans divers regimes. 

5.1 Fonction de partition, energie libre et moments 

La fonction de partition et l'energie libre sont des observables fondamentales. Elles ont ete definies 
precedement pour les trois types de modeles de matrices. 

II est conjecture 1 que la fonction de partition du modele normal est une fonction r isomon- 
odromique au sens de Jimbo-Miwa-Ueno jlHl EH1 E3 EE] , et il est connu que la fonction de partition 
du modele formel est une fonction tau KP [21] . Dans les trois modeles, l'energie libre a une limite 
N grand, reliee a la hierarchie de Toda "sans dispersion" [2*T| I78j. 

Les derivees de la fonction de partition par rapport aux coefficients des potentiels donnent 
acces aux moments: 

tr M l 2 s ) 

9 9kr d~ 9h . . . d- gis Z 

Souvent, il est preferable de calculer les cumulants: 

^con„. ,;,.,,,:/,,,/. := N r+S ~ 2 < tr M* . . . tr Mt tr M l J . . . tr M l 2 s ) conn 

s 

- (-ty +s - 2 HhHk d 9ki . . . d 9kr d~ 9h . . . d~ gis f 

i=\T 1=1 

- 2) 

Dans le modele Normal et Normal-brise, les T conn k , kr .j-. .-^ sont des integrales convergentes. 
Dans le modele formel, les T conn k , kr .j-. sont des fonctions generatrices de graphes ouverts 
j2UE0], avec r bords de signe +, et s bords de signe — . 

5.2 Valeurs moyennes de traces mixtes 

Notons que les derivees de la fonction de partition et de l'energie libre donnent acces seulement 
aux moments non mixtes, i.e. aux valeurs moyennes de produits de traces, telles que chaque trace 
ne contienne qu'un seul type de matrice (Mi ou M 2 ). 
1 Cela a ete prouve dans le cas d'une seule matrice [P7]. 



T, 



N' 



(i i.5 - r 



tr Mf 1 . . . tr M{ r tr M^ 1 . . . 

t ^ rii=i'- rii=i h 

' N 2 



d, 



9k-, 
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Pour les trois modeles de matrices, on peut souhaiter calculer la valeur moyenne de toute traces 
de la forme: 




:= N-EiU / Y[ Tr (m^M 1 ^ 1 . . . M^ ra M l 2 s ' 




(II.5-3) 



dont l'exemple le plus simple est: 



1 



(Tr M*M l 2 ) 



(II.5-4) 



T, 



k,l ■ — 



Ces observables ne peuvent pas etre obtenues en derivant l'energie libre par rapport aux coefficients 
des potentiels. Elles ne peuvent pas non plus etre obtenues a partir des densites de valeurs propres 
(en effet, on ne peut pas exprimer tr M\M l 2 a partir des valeurs propres de M\ et M 2 ). Leur calcul 
est plus difficile, et represente l'un des defis de ce modele. 

Pour le modele normal, l'observable eq. (|II.5-4|) a ete calculee par (Bertola,Eynard) dans [P4], 
et une generalisation a toutes les autres observables mixtes a ete obtenue par (Eynard, Prats- 
Ferrer) dans j3H], en utilisant le theoreme . II 

Dans le modele formel, les observables mixtes du type eq. (|H.5-3|) sont les fonctions generatrices 
de surfaces ouvertes, dont les bords peuvent porter plusieurs signes, autrement dit, avec des 
conditions de bords non-triviales [21]. L'etude des operateurs de bords est un sujet actuel et actif 
en physique des theories conformes et en theorie des cordes, et oil beaucoup reste a comprendre 
[53] . Les modeles de matrices pourraient apporter des resultats a ce domaine. 

Une large classe de ces observables (incluant eq. ([H.5-4j) ) a ete calculee dans la limite N — ► oo 
par (Eynard) dans [P3] et dans [32], et le calcul de toutes les autres observables mixtes est en 
voie d'etre acheve a ce jour. 

5.3 Densites et correlations de valeurs propres 

Considerons ici uniquement le modele normal. Nous allons suivre les traces de Mehta [HHl EH] • 

Donnons nous r et s deux entiers compris entre et N (et l'un au moins non nul), et des 
entiers k%, . . . , k r compris entre 1 et di et l\, . . . , l s compris entre 1 et o^- 



Definition 5.1 Densites de probabilites integrees de eq. \II.2-T^I . sur [Y[i=i^ x ki x n^=i F y i 3 



Pi 



r\s 



(x u . . . 



x r ; yi,..., y s ) dxi... dx r dy x . . . dy. 



s 




dj^Novm, ki,—,kN;h,—,lN 

(xi, . . . ,x N ;y 1} . . .,y N ) 



(II. 5 - 5) 
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C'est la probabilite que X\ G F x kn ■ ■ ■ ,x r E T x k r soient valeurs propres de M\ et y% G T y i 1 , ■ ■ ■ ,y s G 
T y i s soient valeurs propres de M2, simultanement. 

Definition 5.2 Fonctions de correlations de valeurs propres 

r?ki,...,k r ;h,... ,l s ( v V - V V\ 

n r . s [ai, . . . , v\ r , r l5 . . . , r s ) 

j r s 

:= N- r - s ( JJ Tr 5 {xM \Xi - M x ) Y[ Tr 5 (yM {Yi - M 2 ) 

\i=l i=l 
N 



■■="— e e n 

fe r+1 ,..,fe Jv z s+1 ,..,« JV i=i • / **er»* i ■ / wer»j 4 

r / N \ s / N \ 

n E^-'i) n e^-») 

(xi, . . .,x N ;y h . . .,y N ) 

(1 1. 5 - 6) 

o?x (resp. 5 <yV,k \y)) est la distribution 5 de Dirac sur le contour T x k (resp. T y k)- 

On peut aussi definir la partie connexe (les cumulants), on la note: 

^conn'.'^s l '"' ,la (^1) • • • > X r ', Ei, ... , E s ) (II. 5- 7) 

Remarque 5.1 Les deux fonctions i? r;s et p r - s sont presque egales. Lorsque tous les Xi sont distincts 
et tous les yi sont distincts on a: 

R k, r , kr ;h,...,i s . . . ; Xr; yi; . . . ; y s ) = jy ^ rljV _ 8 , /gfe- ,fcr;tl '"" t '(gi. • • • » xr) yi, ■ ■ ■ , y s ) (H.5-8) 

Remarque 5.2 On peut retrouver les quantites definies en eq. (|II.5-l|) : 



T — — 

■ J -pi;...;p r ;qi;—;q a 



E E 



~P1 ^Pr^M „3s 
x l • • • x r ill • • • Us 



Rk},-,kr;h,-,is( Xl: ..., Xr ;y 1 ,...,y s )dx 1 ... dx r d yi ...dy s (II.5-9) 



Remarque 5.3 On a la relation de recurrence evidente: 

Pr-i-s kr ~ 1,lu '"' la (x 1 , x r -D yi,...,y s ) dxi... dx r _xdyi ...dy s 

E / Prff kr '' h ''"' ls (xi, . . . ,x r ;yi, . . . ,y s ) dx\ . ..dx r dyi ...dy s 



k r J xrer* kr 

(II. 5 - 10) 
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5.4 Fonctions generatrices 

II est plus commode d'introduire des fonctions generatrices pour les moments 
Definition 5.3 La resolvante de la matrice Mi est: 

<*<*> == E ^ 4 < tr J±m) = E A * (H6-H) 

fc=0 x 7 k Jy k 

Notons que W\{x) n'est qu'une serie formelle en 1/x (pour les trois modeles), et il faut faire 
attention si on veut l'interpreter comme une serie convergente, et comme une fonction analytique 
de x. En fait W\{x) est une serie asymptotique au voisinage de x — > oo qui definit plusieurs 
fonctions analytiques differentes. 
De meme, 

Definition 5.4 La resolvante de la matrice M 2 est: 

^(z/):=f;^r = ^(tr^— ) = J2 I —,Mf)<V ( IL5 " 12 ) 
yi+i N \ y-M 2 / *Y J TVl y-y' 

Comme precedement, c'est une serie formelle, qui correspond a plusieurs fonctions analytiques 
differentes: 

On introduit aussi des fonctions a deux points (voir notations de eq. (|11.5-2p ): 
Definition 5.5 Fonctions a deux points 

OO 00 rp I \ 

W ll2 (x; y) :=J2E 5^ = ( tr tr ^rjf) ~ ^(,)WM (115-13) 



fc=0 1=0 

oo oo 



fc=0 z=o 

W 2 . 2 (y; y 7 ) = T^'iL = ( tr— — tr —) - N 2 W 2 (y)W 2 (y') (II.5-15) 

t^ ^y k+1 y' l+1 \ y~ M 2 y'-M 2 / 2yy> 2yy ' 1 ; 

Notons ce sont des series formelles en x et y, i.e. des series asymptotiques qui definisent plusieurs 
fonctions analytiques differentes. 

On introduit aussi des fonctions a deux points contenant les traces mixtes du type eq. 1)11.5-4)1 : 



Definition 5.6 Fonction de correlation a deux points mixte 

x k+iyi+i n \ 1 x - M 1 y - M 2 



OO OO rp 1 / 1 1 

W»(x, S ) := E E = Tf ( tr — CT— 1Z ) O 1 ' 5 " 16 ) 



fc=0 Z=0 
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Conclusion 

L'objet de l'etude des modeles de matrices est essentiellement de calculer toutes ces fonctions, en 
particulier dans la limite N — > oo. 
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Chapitre 3 



Methode des polynomes biorthogonaux 
et integrabilite 



Cette methode ne s'applique qu'au modele normal, on peut aussi lui donner un sens pour le 
modele formel, mais nous ne le ferons pas ici. Nous nous placerons desormais dans le cadre du 
modele normal defini dans le paragraphe II I2.3I II y a de nombreuses references sur les polynomes 
orthogonaux: [2B1 E7J HI] , voir aussi [T| 15T] . 

1 Polynomes biorthogonaux 

Definition 1.1 Considerons deux families de polynomes moniques it n {x) = x n + . . . et a n (y) = 
y n + . . ., telles que: 

r dxdyn n (x)a m (y) e -^ v ^ +v ^'^ = h n 5 nm (III.l-l) 



1.1 Generalisation de la formule de Heine 

Theoreme 1.1 Formule de Heine generalisee pour les polynomes biorthogonaux 

J Hn(F) dM 1 dM 2 det{x - M 1 )e-^ M ^ +V ^- M ^ 



= TTT^Z W^TTZTT^TTZ-, 7T7Z\ ( IIL1 " 2 ) 

»CT) 



J H (r) d M 1 dM 2 e-Tiy^)+y^)-^] 



autrement dit, n n est la valeur moyenne du polynome caracteristique de la matrice M\. De meme: 

f Hn(r) dM 1 dM 2 det(y - M 2 ) e~% W^+v^-m^] 

a n(y) = - — — _JV ri/ , u >, IA ,,, , (III. 1-3) 



L m , dM 1 dM 2 e-^ v ^ +v ^- M ^ 



Ce theoreme semble etre du a Jean Zinn-justin [77]. II est expose dans [Pll]. Ecrit avec la 
mesure induite sur les valeurs propres, il donne une formule semblable a celle de Heine jHZ] pour 
les polynomes orthogonaux. 
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II est clair que l'existence des polynomes biorthogonaux est equivalente a la non-annulation 
du denominateur. Pour le modele hermitien, avec des potentiels V\ et V% reels, l'integrand du 
denominateur est strictement positif, et le denominateur ne s'annule jamais (voir aussi |2E])- Pour 
le cas general, il n'a pas encore ete trouve de critere simple sur V\ et V2 pour garantir l'existence 
des polynomes orthogonaux. 

Remarquons d'apres [P6], que le denominateur est un polynome de degre < n dans les k, et 
done la non-existence des polynomes biorthogonaux (i.e. 3n tel que le denominateur s'annule) est 
un sous-ensemble de mesure nulle (denombrable) de C dlXd2 . Autrement dit, pour V\ et V2 donnes, 
les polynomes biorthogonaux existent pour presque tout k. Cet argument est du a M. Bertola. 

1.2 Notation, fonctions d'ondes 

Nous allons suivre les notations de [P8], [P6]: 
Definition 1.2 Fonctions d'ondes 

V "n V "n 

Definition 1.3 Transformees de Fourier- Laplace des fonctions d'ondes 

k = l,...,d 1 , fl*\y):= [ Mx)^ Xy dx 

k = l,...,d 2 , [ Mv)^ xy dy (HI.1-5) 

Definition 1.4 Transformees de Hilbert des fonctions d'ondes 



V>W(i/):=e^w / ^l e -^>e^ xy dxdy' (III.1-6) 

Jvv-y 



r x — x 



(III.1-7) 



Notons que les tp n , <p n , ainsi que les ip^ et (p^ pour k > sont des fonctions entieres. Par 
contre, les transformees de Hilbert sont analytiques par morceaux dans C/UjT x j (resp. C/UjT y j). 
La fonction <f>n\x) est discontinue chaque fois que x croise un contour T x j, et sa discontinuity 
vaut (x + et x_ designent respectivement des points a gauche et a droite du contour): 



5f ) (x + ) = #?)(x. 



+ 2iirJ2 K i,i<t>n( x ) (III. 1-8) 



1=1 

de meme, si y £ r^,-: 

A 0) (y+) = ^ ) (y-) + 2i7rJ2^,j^ ) (y) ( IIL1 -9) 
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On peut aussi definir des fonctions entieres a partir des transformers de Hilbert, par prolongement 
analytique en deplagant les contours, pour chaque composante connexe A de C/U,- T x j, on obtient 
une fonction entiere sur C, qui coincide avec ipn\x) sur A. Cette procedure permet d'obtenir d\ 
(resp. d 2 ) fonctions entieres independantes, notees: 

$$*\x) , k = 1, . . . , di , ^> k \y) , k = 1, . . . , d 2 (III.1-10) 



2 Relations de recurrences et matrices a bandes finies 

Les polynomes 7r n forment une base, et Ton peut decomposer le polynome xn n (x) sur cette base, 
et done pour les ip n \ 

n+l 

xtp n (x) = ^2Q nm ip m (x) (III.2-1) 

m=0 

de meme: 

n+l 

y<t>n{y) = Y. p ™^y) (in.2-2) 

m=0 

ou Q et P sont des matrices semi-infinies. 

De meme, on peut decomposer le polynome n' n (x) sur la base des n n , le terme de plus haut 
degre etant n7r n _i. La derivee de ip n peut done se decomposer sur la base des tp m , avec m < n + dx 
(a cause de la derivee de e~^ Vl ). Les coefficients de ce developpement sont obtenus par l'identite 

Vn,m, J^dxdyd x {f n {x)<f> m {y) e^) = (III.2-3) 

qui implique: 

rp n+di 

--d x i) n { X )= P mn^m{x) (111.2-4) 

m=n— 1 

On a vu plus haut que cette somme est en fait limitee a m < n + d±, i.e. 

Pnm = si m>n+loum<ti - deg V[ (111.2-5) 

de meme: 

Qnm = si m>n+loum<H - deg V 2 ' (III. 2-6) 

i.e. les matrices Q et P sont a bande finie. 
On pose: 

ttfe(n) := Q n ,n-k , /3fc(n) := P n ,n-k (III.2-7) 
Le fait que les polynomes sont moniques implique: 



n_ 1 (//)= X l^:- ~„ (TTI.2-8) 
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Autrement dit: 



Q 



( «o(0) 
«i(l) 

a d2 (d 2 ) 


V o 

/ A)(o) 



7o . . . 

«o(l) 7i 



a d2 (rf 2 + 1) 






ld 2 



7o 








(III.2-9) 



\ 



A)(l) 7i 



/3 dl (di + l) 

; o 
V o 





7di 







(III.2-10) 



2.1 Notations pour les matrices semi— infinies 

Si on definit la matrice semi-inflnie: 



A := 




i.e. A,-,,- = Si 



ij — 



(III.2-11) 



et les matrices: 



a k := diag(a fc (n)) n=0> .. 
on peut ecrire: 



f3 k := diag(/3 fe (n)) n=0r ..,oo , 7 := diag(7 n ) n=0> ... )OO (III.2-12) 



d 2 



di 



Q = 7 A + ^a fc A* fc , P = 7 A + ^^A' fc (III.2-13) 

k=0 k=0 

On introduit aussi les projecteurs sur le sous espace des polynomes de degre < n: 



On a les relations: 



Il n := diag(l, . . . , l(n + 1 termes), ...,1,0,...) 



An n = n„_!A , aa* = i , a'a = i - n 



(III.2-14) 



(III.2-15) 



Remarque 2.1 Notons que l'algebre des matrices infinies considered ici, n'est pas associative (A(BC) ^ 
(AB)C), car on ne peut pas echanger l'ordre de sommation pour les sommes infinies. II n'y a associativite 
que si toutes les sommes impliquees sont finies, c'est le cas pour les matrices a bandes finies et les 
projecteurs H n . 
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2.2 Relations entre P et Q 

Le fait que la decomposition de n' n sur la base des n commence par: 

ir' n (x) = nir n -i(x) + . . . (III.2-16) 

implique pour P et Q les relations: 

Tn 

(P%m = (Vl(Q)) n , m Si m > 71 , (P*)^-! = (K(Q))n,n-l - T7 (III.2-17) 

iV7n-l 

De meme: 

(Q*)n,m = (V 2 '(P))n,m si m > U , (Q')n, n -1 = (^Wkn-l " "TT (III.2-18) 

iV7n-l 

Ces relations sont souvent appelees "equations du mouvement" [2*Tj . 
Cela implique en particulier: 

A d *a d2 = ~g d2+1 (iA) d2 , A dl P dl = g dl+ i(lA) dl (III.2-19) 
On a aussi la relation de Heisenberg aussi appelee "equation des cordes" [2*T] : 

[P\Q} = ^1 (III.2-20) 

2.3 Relations de recurrences, resume 

On a l'ensemble de relations: 



(III.2-21) 



-jfd x ip n = Yfk=-iPk{n + k)i) n+k , -jfd y (j) n = Y?k=-i a ki n + k)(j) n +k 
De plus, il est facile de montrer [P6] que pour tout j G [0, d 2 ]: 

\ / 0,n 

et pour tout j e [0, di]: 



iV 



aw? = Et 1 « i Wa + ^*l'-(.)(«H) < nl ' 2 - 23 > 
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3 Noyaux et densites 



On introduit [P10]: 
Definition 3.1 les noyaux 



N-l 



K 12 (x,y):=J2M x )<f>j(y) ( m -3-l) 
i=o 

N-l 

K^{x,x') ■^Y.^W^') ( IIL3 - 2 ) 

3=0 
N-l 

kSV.v) ■■=J2^ k) (y') ( f>M ( IIL3 - 3 ) 

j=0 



di d 2 N-l 

K 21 ( y >,x') : = £X>><E^fV 

k=l 1=1 3=0 

di d 2 oo 

= -EX><E^fV)0?V) (in.3-4) 

k=l 1=1 j=N 

Ces noyaux permettent de calculer les densites definies en eq. (|II.5-5|) . 

Theoreme 3.1 (Eynard, Mehta [P10]) les densites de valeurs propres definies en def.H T5~l\ sont 
donnees par: 

'(.r,. . . . , x r ; Vl , . . . , y s ) = det ( *$) Xi > ^ V \ ) (III.3-5) 

V K 2 i{yi,x j ) K^ 2 \y h yj)J 

La preuve est exposee dans [P10], c'est une generalisation du theoreme de Dyson pour les 
polynomes orthogonaux Elle utilise la relation de recurence eq. fjll.5-10|) et les proprietes 

auto-reproduisantes des noyaux: 



Kk,k> 

k,k> J 


j x du K^\x,u] 
Tk 


)*$?« 


[u,x r ] 


) = KS(x,x') 


(III.3-6) 


Kk,k> 
k,k' J 


/ duK { 22 \y,u 


)K$ 


(u,y' 


) = *&W) 


(III.3-7) 


k,l 


L du K 21 (y,u) 
JVk 


K$( 


u, x) 


= K 21 (y,x) 


(III.3-8) 


k,l 


\ duK [ 22 \y,u 


)K 21 { 


u,x) 


= K 21 (y,x) 


(III.3-9) 



Grace a ce theoreme, le calcul de ces quatre noyaux (notons que ceux ci s'obtiennent par trans- 
formee de Fourrier a partir de Ki 2 ), est suffisant pour calculer toutes les fonctions de correlations 
de valeurs propres du type eq. f)II.5-5|) . 



31 



4 Matrices de Christoffel— Darboux 

Theoreme 4.1 La matrice: 

A n :^[Q,U n _ 1 ] (III.4-1) 
possede settlement un sous bloc de taille d 2 + 1 x d 2 + 1 non identiquement nul: 

{&k-i(k) si n<k<n + d 2 — 1 , n — d 2 <l <n — 1 , k — I < d 2 
-7„_i si /c = n- let/ = n (III.4-2) 
autrement 

De meme, B n := [P, n n _i], possede seulement un bloc de taille di + 1 non nw/: 

(Pk-i(k) si n<k<n + d\ — 1 , n — d\ < I < n — 1 , k — I < d\ 
-7„_i si fc = n- letZ = n (III.4-3) 
autrement 

Definition 4.1 Les matrices A n et B n seront appelees matrices de Christoffel— Darboux. 

On a immediatement la generalisation du theoreme de Christoffel-Darboux [P8]: 
Theoreme 4.2 theoreme de Christoffel-Darboux generalise 

(x'-x)K^ n (x,x') = E £ A n ^WfV) (III.4-4) 

i=n— 1 j=n—d,2 
n+di— 1 n 

(y'-y)<n(y',y) = E E BrnJPMMy) ( IIL4 -5) 

i=n— 1 j=n—d\ 

Grace a ce thoreme, le calcul des noyaux se reduit au calcul de seulement d 2 + 1 polynomes de 
type ■?/>„, et di + 1 polynomes de type <f> n . Ceci est particulierement utile dans la limite n — > oo. 
II n'est pas clair qui a decouvert ce theoreme pour la premiere fois. 



Corollaire 4.1 On a des relations additionnelles (Bertola, Eynard, Hamad [P8]j, par exemple: 

n+d2 — 1 n 

)M\ . (x. ii \ - 



(x + ^d y )K 12 , n (x,y) = - E E A ni ^ 3 (x)Uy) (HI.4-6) 

i=n— 1 j=n—d'z 



rp n+di — 1 rt 

(y + ^0 x )tfi2,„(*,y) = - E E B T*j1>i(*)Mv) ( IIL4 " 7 ) 

i=n—l j=n—di 
rp n+di — 1 71 

^(4 + ^)^!jx,x')= E E a*jiM*)# fc V) (ni-4-8) 

i=n—l j=n—di 
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5 Repliement sur une fenetre 



On a vu que les noyaux peuvent s'exprimer seulement a l'aide des polynomes de rang £ [n — d 2 , n] 
et des transformees de Fourrier de rang e [n — l,n + d 2 — 1]. 

Definition 5.1 Les vecteurs suivants, formes de d 2 + 1 fonctions d'ondes ip (resp. transformees 
de Fourrier <j)) consecutifs sont appelles fenetre et fenetre duale: 

¥„(*) := (iP n „ d2 (x),...,iP n (x)Y , *$\x) ■= (0j£i(^),---^J2*-i(*))' (IIL5-1) 
et de meme: 

Definition 5.2 Les vecteurs suivants, formes de d\ + 1 fonctions d'ondes <fi (resp. transformees 
de Fourrier ip) consecutifs sont appelles fenetre et fenetre duale: 

*»(V) := (<j>n- d M,---,Uy)y , ¥ n k \y) := (^(y), . . . , ^^(y))* (IIL5-2) 

II est clair qu'on peut utiliser la premiere relation eq. (jIII.2-21j) (multiplication par x) pour exprimer 
n'importequel polynome ip m avec m e [0, oo[, comme une combinaison lineaire a coefficients 
polynomiaux en x de polynomes de la fenetre. 

Definition 5.3 La matrice F n (x) de taille oo x (d 2 + 1) telle que 

n 

Vm = 0,...,OO, ^ m {x)= Yl ( F n( X ))m,jM X ) ( IIL5 " 3 ) 

j=n— d.2 

est appelee matrice de repliement sur la fenetre \l/ n . 

On a deg {F n ) m - < n — d 2 — m si m < n — d 2 et deg {F n ) m - < m — n si m > n. La matrice F n (x) 
a ete calculee par Bertola-Eynard dans [P2], et vaut: 

Theoreme 5.1 (Bertola, Eynard [P2]J La matrice F n (x) est donnee par 

F n (x) = ((Q - x)^ - (Q - x)^ 1 ) A n (III.5-4) 

oil (Q — x)^ 1 et (Q — x)^ 1 sont respectivement les inverses a droite (triangulaire inferieure) et 
a gauche (triangulaire superieure) de la matrice Q — x. On les calcule en introduisant les deux 
matrices suivantes respectivement strictement triangulaires superieures et inferieures: 

Q L :=l-J— (7 _1 A) d2 (Q-x), Q R := 1 - n - AV 1 (Q - x) (III.5-5) 

9<h+i 

on a: 

(Q - ^z 1 = ^- (i - Ql)- 1 (7 _1 A) d2 (HI.5-6) 
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(g-x) i? 1 = (l-Q /? )- 1 AV 1 (HI.5-7) 

la matrice (1 — Ql)~ 1 est une notation pour YlT=o(^L) k '■> e ^ ^ ^ au ^ n °ter 1 ue pour le calcul d'un 
element de F n (x), cette somme se reduit a une somme finie (cf remarque III 12.1|) . 

De meme on definit les matrices respectivement triangulaires inferieure et superieure: 

Q L :=1- n d2 _! - J— (AVY 2 (Q - x)* , Q R :=1- 7 _1 A (Q - xf (III.5-8) 

9d 2 +i 

et la matrice: 

F n (x) = (J- (1 - Ql)- 1 (AV 1 )* - (1 - e«)- 1 7" 1 A) 4 (III.5-9) 

5^2 + 1 

qui assure le repliement des transformers de Fourrier (pour tout k = 1, . . . , d?): 

n+d 2 -l 

Vm = 0, . . . , oo , £ F mj>n (x) ~^f\x) (III.5-10) 

j=n-l 

Et on peut definir de fagon similaire les repliements des <fi et ip: 

n n+d\ — 1 

<f> m (y) = E GnnuMMv), 1>ZHv)= E ^(y)^ fc) (y) (ra.s-ii) 

j=n-d 1 j=n-l 

6 Systemes different iels 

Theoreme 6.1 La fenetre ^ n (x) satisfait un systeme d'equation differentielles lineaires a coeffi- 
cients polynomiaux 

- ^d x * n (x) = V 1>n (x)* n (x) (III.6-1) 
est /e sen/ non nw/ la matrice: 

v x , n {x) = (n„ - n n _, 2 _!) p< (m.6-2) 

preuve: (voir [P8], [P2], mais ces systemes ont ete trouves dans des cas particuliers avant) 
on utilise le repliement sur une fenetre. On ecrit la relation eq. (jIII.2-21|) (d x ip) pour tous les 
polynomes de la fenetre, et Ton utilise eq. (|III.5-3|) pour exprimer le membre de droite sur la base 
des polynomes de la fenetre: 

rp m+di m+di n 

Vm = 71 - d 2 , . . . ,71 , - — d x 1p m (x) = E Pkm^k{x) = E E P kmFkj,n{ X )i ) j( X ) 

k=m—l k=m—lj=n—d2 
n 

= E (P'Fn^mjlpAx) (HI.6-3) 
j=n— d2 

QED. 
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Notons que 2\ n (x) est une matrice de taille 4 + 1x4 + 1, dont les coefficients sont des 
polynomes en x de degre < d\. 

De meme, on definit trois autres systemes differentiels (on suppose n > max(c?i, d 2 ), sinon il y 
a des termes correspondant aux petites valeurs de n comme dans eq. (jlll.2-22j0 : 



Definition 6.1 Autres systemes d 'equations differentielles lineaires a coefficients polynomiaux 

v hn (x) = (n n+d2 _x - n n _ 2 ) PF n (x) (in.6-4) 



^dM\x)=V 1 ,n(xM'^> 



T 



- jjdy<$> n {y) = V 2>n {y)$ n {y) , T> 2 , n (y) = (U n - n n _ dl _x) Q l G n (y) 



T 

N 



d y ¥ n k \y)=V 2tn (y)¥ n k \y) , V 2>n (y) = (U n+dl ^ - Il n _ 2 ) Q G n {y) 



(III.6-5) 
(III.6-6) 



Notons que t>i iTl (x) est une matrice de taille d 2 -\-lxd 2 + l, dont les coefficients sont des polynomes 
en x de degre < di, comme P l n (a;). De meme, T> 2n {y) et X\ n (2/) sont des matrices de taille 
d\ + 1 X di + 1, dont les coefficients sont des polynomes en y de degre < d 2 . 

Dans [P2], nous avons obtenu des expressions encore plus explicites de ces systemes. 

Theoreme 6.2 (Bertola, Eynard [P2]j La matrice T>\ )n {x) s'ecrit: 

• Vr(Q)n-fc,n-l \ 



( V{(Q)n-d 2 ,n-d 2 ■ 

o '■. ; 

V/(Q)n-l,n-l 





+ 



'7(71 - d 2 - 1) 



/ V{(Q)-V{(x) 

Q~ x n—d,2,n—l 

V{{Q)-V{{x) 
\ Q~ x n,n—l 







7(n - 1), 



/ 



V 



o 

Qd 2 -l("-l) 
f>d 2 ("-1) 







x— ao(n— 1) 



7("~!) \ 
ad 2 («-l) I 






V{{Q)-V{{x) ^ 
Q—x n —d 2 ,n+d2 — l 



V{(Q)-V{(x) 



A N 



[111.6-7) 



n,n+d2 — l 



J 



La preuve est detaillee dans [P2]. Elle utilise ea . fllll.6-21) et ea. 1111.2-171) . 
On a bien sur des expressions similaires pour les 3 autres systemes [P2]. 



7 Dualite spectrale 



II a ete pouve par Bertola-Eynard-Harnad [P8] que ces systemes differentiels sont duaux. 
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Theoreme 7.1 (Bertola, Eynard, Hamad \P&\) les matrices 2\ n (x) etT>{ n (x) sont conjuguees 
par la matrice de Christoffel-Darboux: 

A n V 1>n (x) = V\ >n (x)A n (III.7-1) 

preuve: donnons ici une idee de la demonstration de [P8]: 

Supposons que n > d 2 + d\. Considerons deux solutions arbitraires des systemes T> ln {x) et 
i>i,n{x)- 

f(x) = (f n . d2 (x), f n (x)Y , ^dj(x) = -V hn (x)f{x) (III.7-2) 

T 

g(x) = (9n-i{x), . . . ,flr T1+<fa _i(z)) , —d x g(x) = Vi >n (x)g(x) (III.7-3) 

On definit les fonctions f m {x) avec n — d 2 — 1 < m < n + d\ et g m (x) avec n — d\ — 1 < m < n + d 2 , 
recursivement par: 

^ otk{m)f m -k{x) = xf m {x) , a fc (m + k)g rn+k { x) xg m (x) (III.7-4) 

fe=-i fc=-i 

les equations eq. (|III.7-2|) et eq. ()III.7-3|) sont alors equivalentes a (voir [P8] pour details): 

rp dl 

Vm e [n - d 2 , n] - —d x f m (x) = ^ Pk(m + k)f m+k (x) 

k=-l 

rp dl 

Vm G [n - 1, n + d 2 - 1] —d x g m {x) = ^ Pk(m)g m ^ k (x) (III.7-5) 

fe=-i 

c'est a dire les relations eg. f)III.2-21[) et eq. (jIH.2-22|) . On a done: 

T 

N 



g'ix') (v[ >n {x')A n - A n V x , n {x)) f(x) 

jf E E ftW.oo (^U(x'Mn-A n p lin (x)) r 

r=n—l s=n—d2 

]T £ g r {x')f s {x)[P\A n ] rs 

r=n—di — l s=n—d'2 — l 
n+cfa n+d\ 

E E 9r(x')f s (x)[P\[Q,U]] r)S 

r=n—di — l s=n—d'2 — l 
n+d2 n+di 

E E 9r(x')f s (x){[[P t ,Q],U] + [Q,[P t M) i 

r—n—di—1 s=n—d2 — l 
n+d2 n+d\ 

E E 9r(x')fs(x)[Q,[P\n]\ r)S 

r=n—di — l s=n—d2 — l 

n n+rfi— 1 

V~ X ) E E 9r{x')f s {x)[P\U] rs 
r=n—d\ s=n— 1 
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n n+di — 1 
(X-X') £ E 9r(x')f s (x){Bi) rs 
r=n—di s=n—l 

(J//.7-6) 

qui s'annule pour x' = x. Ceci etant vrai pour toute solutions / et g, on doit avoir: 

Vl n (x)A n = A n Vi, n (x) (III.7-7) 
QED. Notons que Ton a utilise la relation de Heisenberg entre P et Q: 

[P\Q] = Li (HL7-8) 

A partir de la, il a ete pouve par Bertola-Eynard-Harnad [P8] que ces quatres systemes 
differentiels ont la meme courbe spectrale: 

Theoreme 7.2 (Bertola, Eynard, Hamad [P8]J Les 4 systemes T>i tn (x), V ln (x), T> 2 , n (y), T> 2>n (y) 
ont la mime courbe spectrale: 

E n (x, y) = g d2+ i det {yl - Vi >n (x)) = g d2 +i det [yl - V^ n {x)\ 

= 9d x +x det (xl - T>2, n (y)) = gdx+x det (xl - T>2, n (y)j (III. 7-9) 

idee de la preuve: les egalites diagonales (par exemple entre 2\ n (x) et T>2, n {y)) ne font pas 
intervenir les proprietes d'orthogonalite, elles sont consequences de la transformee de Fourrier, et 
peuvent etre obtenues par simple reecriture des deux membres [P8]. 
Les egalites horizontales sont consequences du theoreme III I7.1I 



8 Solution fondamentale 

Considerons ici le systeme T> 1>n (x), et supposons n > d\. Nous savons deja que chacune des 
c?2 fonctions $n\x) avec k = 0, . . . , d 2 est solution de ce systeme. La solution fondamentale de 
T>i n (x) a ete trouvee par Bertola-Eynard-Harnad dans [P6], et simultanement par A. Kapaev 
•")() seulement pour le cas d\ = d 2 = 2. 

Theoreme 8.1 (Bertola, Eynard, Hamad [P6]J La matrice carree *& n ( x ) dont les colonnes sont 
les vecteurs <f>J?(x) avec k = 0, . . . , d 2 , est une solution fondamentale du systeme T> l n (x): 

^ d x * n (x) = T> 1>n (x) # n (x) (III.8-1) 

preuve: Pour n > di, les relations de recurrences eq. (jIII.2-22|) sont les memes pour j = que pour 
j = l,...,d 2 , ce qui implique que la fenetre duale &n\x) (ou <fin\x) designe l'une quelquonque 
des fonctions entieres defmies par eq. (fill. 1-6(1 ^1 est aussi solution du systeme P l n (x), et done $„(x) 
satisfait eq. (jlll.8-l|) . Pour verifier que cette matrice est bien une solution fondamentale, il suffit de 
calculer son determinant, i.e. le Wronskien, et verifier qu'il ne s'annule pas. Le theoreme 111 16.21 
implique (voir [P2]): 
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Lemme 8.1 Trace de la matrice T>i %n (x): 

tr V ltn (x) = V((x) - J^- (III.8-2) 

9d 2 +i 

ce qui implique 

^- In det $ n (x) = tr ^(xf 1 ^-^) = tr V 1>n (x) = V((x) - J*- (III.8-3) 

dx g d2+1 

et done 

N_( -tr / \_ ^d 2 \ 

det $ n (x) = C e T V Sj s+ 1 / (III.8-4) 

La constante C a ete calculee dans la limite x — > oo dans [P6], [P2], et elle a ete calculee 
directement dans [P2]. Elle est non nulle. QED. 

Avant d'etudier les asymptotiques a x grand dans le paragraphe suivant, mentionons les solu- 
tions fondamentales des trois autres systemes. La solution du systeme 2\ n (x) a ete explicitee par 
[P6] sous une forme recursive, et explicitee precisement pour certains cas particuliers: potentiels 
cubiques par [SD], et un cas particulier de potentiels de degres d\ = 7, di = 11 par [P6]. II apparait 
clair que la construction de [P6] est generale, mais une expression fermee reste a trouver. 

On peut aussi construire une solution fondamentale de T>i >n (x) en utilisant la dualite 
ea. lllll.7-lj) . 

Theoreme 8.2 (Bertola, Eynard, Hamad [P6]J La matrice 

* n (x) := A- 1 fa^y 1 C (III.8-5) 

ou C est n'importe quelle matrice constante inversible independante de x, est solution de 
-^d x ^ n {x) = Vi >n {x)V n (x). 

En effet, si *„(x) est solution du systeme V 1>n (x), alors on a, d'apres le theoreme III IT. II 

— (** (z)4,¥ B (z)) = -& n {x){V\ in {x)A n - A n V hn (x))* n (x) = (III.8-6) 

QED. 

9 Asymptotiques x — > oo et phenomene de Stokes 

Les asymptotiques de la solution fondamentale & n (x) lorsque x — ► oo ont ete trouvees par Bertola- 
Eynard-Harnad [P6]. 
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Definition 9.1 On introduit les notations suivantes: 

yi(x), . . . , yd 2 {x) son t les d 2 solutions de I 'equation V 2 (y) = x 



T k (x) ■.= J2-^Tj9i+iyk(x) l+1 , T(x) = di&g(V 1 (x),T 1 (x),...,T d2 (x)) (III.9-1) 



uj:=e^, ^oo = l, ttij = uj ij sii,j>0 (III.9-2) 
^ , 2n + 1 — do 2n — 1 + d 2 . 

G = diag(-n, 2d^' " ' ' 2dT~^ (IIL9 " 3) 

n + 1 + 

F = diag(v^, . . . , ) (III.9-4) 

Theoreme 9.1 (Bertola, Eynard, Hamad [P6]J Asymptotiquement pour x — > 00, on a: 

* n (ar) ~r (l + O(a;~ 1/d2 ))a; G fie^ T{a;) C(x) (III.9-5) 

cm C(x) es£ ime matrice constante par morceaux donnee dans [P6]. C(ar) es£ constante dans des 
secteurs d 'angle n /d 2 . 

preuve: Nous n'allons que survoler la preuve detaillee de [P6] . 

Rappelons que par leur definition Hi ll. 21 les fonctions d'ondes 4>j(y) sont entieres, et leur 
comportement asymptotique a grand y est donne par: 

fay) ~ ie-? % W(l + 0(l/y)) (III.9-6) 



Nous calculons alors les asymptotiques a grand x des 0J- (a;) par la methode du col. 
Pour x donne, posons: 

/x(2/):=V 2 (y)-^ (III.9-7) 
On cherche les points cols, solutions de d y f x (y) = 0, i.e. 

V 2 \y) = x (III.9-8) 

II y a G? 2 points cols 3/1(2;), ... , yd 2 (x). Pour = 1, . . . , d 2 , on definit le chemin de col l v k(x): 

"k(x) = f-'ifMx)) + R+) (III.9-9) 

Les d 2 chemins ainsi definis forment une base homologiquement equivalente a la base des Y y k- H 
existe done une matrice de changement de base M(x): 

d 2 

f ; = ^M ifc (x)7 fc (x) (III.9-10) 
fc=i 
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La matrice M(x) a coefficients entiers est constante par morceaux sur des domaines du plan 
complexe. Cette matrice a ete calculee explicitement dans [P6]. Pour chaque chemin 7\(o;), on 
a 1' approximation du col: 

/ My)e^dy ~ V -M, J NV %, e^^l + 0{x~^)) (III.9-11) 

et done: 

~ E M, t ^g ^ .-**<*(■»(! + 0(x-^)) (IH.9-12) 

En ce qui concerne l'asymptotique de ^\x), choisissons x — ► oo dans un secteur qui ne croise 
pas de T x k- Dans ce secteur, la relation d'orthogonalite implique: 

4>f\x) ~ y^x-^+V e ^ x \l + O(x^)) (III.9-13) 

Les asymptotiques dans les autres secteurs se calculent par la relation de discontinuity eq. (fill. l-8j) . 

La matrice C(x) est en partie constitute de la matrice M(x) definie en eq. (|III.9-10|) . et des 
relations de discontinuite eq. (|III.l-8j) . Elle est donnee precisement dans [P6]. QED. 

Les disco ntinuites de la matrice C(x) representent les matrices de Stokes: 

C(x^)- 1 C(x+) (III.9-14) 

Celles ci sont calculees dans [P6], leur expression prendrait trop de place ici. Ce qui est important, 
e'est de remarquer que les matrices de Stokes, de meme que les conditions de saut de eq. (|III.l-8|) 
sont independantes de n, de x, et des coefficients des potentiels. Nous avons done affaire a un 
probleme isomonodromique. 

10 Probleme de Riemann— Hilbert 

Posons nous le probleme suivant: 

• existe-t-il une matrice *& n (x) analytique par morceaux, et inversible pour tout x, qui satisfait 
les asymptotiques eq. fill. .9-15)) . et les conditions de saut eq. MII.l^l ? 

C'est le probleme de Riemann-Hilbert. La reponse est oui est la solution est unique. 

Le probleme de Riemann-Hilbert pour les polynomes biorthogonaux a ete trouve simul- 
tanement par A. Kapaev [SU] pour le cas ou d\ = d 2 = 2, et par Bertola-Eynard-Harnad [P6] 
pour le cas general. Notons qu'un autre probleme de Riemann-Hilbert a ete etudie dans [2*%] . 
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L'approche par le probleme de Riemann-Hilbert est particulierement utile si Fori s'interesse a 
la limite n grand. En effet, les discontinuity eq. (|III.l-8jl sont independantes de n, et les asymp- 
totiques eq. (jIII.9-5|) en dependent d'une fagon tres simple. 

L'idee est la suivante. On introduit un ansatz ^ n ,ansatz(^) ; dont on connait toutes les discon- 
tinuity et asymptotiques. Si Ton peut montrer que pour x —* oo: 

^n.ansatz^)^; 1 ^) = 1 + 0(l/n) (III.lO-l) 

et que les discontinuity des asymptotiques de $ njansa t z (x)$~ 1 (x) sont d'ordre 0(l/n) partout, 
alors on a uniformement pour tout x: 

K^Ux)^ 1 ^) = 1 + 0(l/n) (III.10-2) 

Cette approche a ete developpee avec un immense succes par Bleher et Its [TU] , et par [TBI CHI H3 
120]. 

Ceci reste a faire pour les polynomes biorthogonaux, mais on peut esperer que ce sera fait tres 
bientot jQE]. 



11 Deformations isomonodromiques et integrabilite 



Remarquons que la propriete d'isomonodromie (i.e. le fait que les matrices de Stokes, et les 
matrices de saut eq. (jlll.l-8|) ne dependent ni de x, ni de n, ni des potentiels implique que la 
matrice ^d x & n (x)\ $ n (x) _1 n'a aucune discontinuite, elle doit done etre analytique entiere, et 
du fait de son comportement a l'mfini, ce doit etre une matrice polynomials Elle doit etre egale 
a Vi )Tl (x) donnee dans eq. ((111.6-7(1 . 

De meme, la propriete d'isomonodromie implique que la matrice Uk, n {x) = 
dg k & n (x)^ $ n (x) _1 n'a aucune discontinuite, elle doit done etre analytique entiere, et 
du fait de son comportement a l'mfini, ce doit etre une matrice polynomiale. De meme pour les 
matrices Vk,n( x ) = \^9g k ^ n (x)j $ n (x) _1 et R n (x) = ^$ n+i (x)j $ n (a;) -1 . Leurs expressions 
explicites sont donnees dans [P2]. 

Remarquons que l'on doit avoir les relations de commutation suivantes: 

T 



N 



T 

N 
T 



-^d x -V 1>n {x),—d 9k 



N 



d x - T>iJx) 



V k ,n{x) 

T 








V^ n+ i(x)R n (x) = R n (x)V 1>n (x) + —d x R n (x) 

T 

Uk,n+i{x)R n (x) = R n (x)U k ,n(x) + —d gk R n (x) 

T 

Vk,n+x{x)Rn{x) = R n {x)V k ,n{x) + —d gk R n (x) 
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(III.11-1) 



Ces relations de compatibilite de Frobenius, signifient que l'on a affaire a des systemes integrables. 

Tout ceci donne lieu a l'existence d'une fonction r au sens de Jimbo-Miwa-Ueno [IHJ|in|, qui est 
conjecturee etre exactement la fonction de partition du modele de matrice Z, definie en eq. (jll.2-14|) . 
comme cela a ete prouve pour le modele a une matrice dans [P7] par Bertola-Eynard-Harnad. 

L'approche integrabilite des modeles de matrices et polynomes orthogonaux a ete abondament 
etudiee [H3 EE E01 113 El • 



12 Traces mixtes 

Nous avons mentione au paragraphe 15.21 qu'un probleme interessant est de calculer des observables 
mixtes du type: 

(trikffil4) (III.12-1) 

A priori, les valeurs propres de M\M\ ne peuvent pas s'exprimer en fonction des valeurs propres 
de Mi et de M 2 , et il semblait improbable de trouver une expression de eq. (jIII.12-l|) en termes 
de polynomes biorthogonaux. Toutefois cela est possible, en utilisant la formule de Morozov du 
theoreme 11 11.31 II a ete trouve par Bertola-Eynard [P4]: 

Theoreme 12.1 (Bertola, Eynard, [P4]J La fonction generatrice de correlations mixtes 
eq. \TTT. est donnee par 

Ji ( tr ^M^) = 1 - de ' ( lN - ^-'^Q^P. 11 "- 1 ) (IIU2 - 2 > 

La generalisation de cette formule pour toutes les autres fonctions de correlations mixtes (i.e. 
produit d'un nombre quelconque de traces, chacune contenant un produit quelconque de matrices 
Mi et M2) a ete trouvee tres recement par Eynard, Prats Ferrer j^Hj, en utilisant le theoreme 
11 11.11 mais sortirait trop du cadre de cette habilitation. 
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Chapitre 4 



Methode des boucles et Geometrie 
algebrique 



La methode des boucles exploite les equations de Schwinger-Dyson, i.e. l'invariance de l'integrale 
par les reparametrisations infinitesimales [2U HSU EHj • Contrairement a la methode des polynomes 
biorthogonaux, celle ci ne fait pas appel a l'integrale d'ltzykson-Zuber-Harish-Chandra. 

La methode des boucles s'applique au modele formel comme au modele normal. C'est toutefois 
pour le modele formel qu'elle est le plus utile, car elle se developpe facilement en puissances de 
1/N 2 . En particulier, la methode des boucles donne facilement la limite N — > oo. 

Dans la limite N — > oo, la methode des boucles fait apparaitre une courbe algebrique, qui n'est 
rien d'autre que la limite de la courbe spectrale des polynomes biorthogonaux vue en eq. (jIII.7-9l) . 
Toutes les observables du modele, dans la limite N — > oo, s'expriment done a partir de proprietes 
geometriques de cette courbe. C'est le lien entre modeles de matrices et geometrie algebrique. 

1 Equations de boucles 

1.1 Changements de variables matriciels 

On considere un changement de variable du type (5 infinitesimal): 



oil / est une fonction analytique de ses deux variables (en pratique, ce sera un polynome). 
La mesure dM\ est multipliee par un Jacobien: 




(IV.1-1) 



dM[ = \J(f(M 1 ,M 2 ))\dM, 




Definition 1.1 La variation Ji(f) du Jacobien a I'ordre 1 en \8\, est definie telle que: 



\J(f(M 1 ,M 2 ))\= (l + 2Re (5 J 1 (f(M u M 2 ))) + 0(\5 2 \)) dM x 



(IV.1-3) 
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De meme, Taction e t tr [Vi(m 1 )+V2(m 2 ) MiM 2 ] eg ^ c h an g£ e) a l'ordre 1 en 5 par: 



e T 



^ tr [V 1 (M[)+V 2 (M 2 )-M[M 2 ] 



(IV.l - 4) 



I 2 Re( ^tr[/(M 1 ,M 2 )(y/(M 1 )-M 2 )] ) +0(|5 2 |) ) e -4 r "' i r ''- W| )+1 - ( - w - ,, ~- Wl - w -'J 



L'invariance de l'integrale par changement de variable Z = J dMidM 2 e T tr [ yi ( Ml )+ y2 ( M2 ) M i M a] _ 
/ dM[dM 2 e-^ tT [Vi(m[)+v 2 (m 2 )-m[m 2 ] implique: 

J dM 1 dM 2 \^J 1 (f(M 1 ,M 2 )) -jtr [/(M 1 ,M 2 )(V 1 , (M 1 ) - M 2 )]j e~ NtI [Vi(a/i)+v 2 (m 2 )-a/ 1 m 2 ] = Q 

(IV. 1-5) 

i.e. 

Theoreme 1.1 equation de boucles: pour tout changement de variables f , on a: 

^ (Ji(f(Mi,M 2 ))) = I (tr [/(M 1 ,M 2 )(^(M 1 ) -M 2 )]> (IV.1-6) 

L'equation eq. (jiV.l-6j) constitue la forme la plus generale possible de ce que Ton appelle les 
equations de boucles. II s'agit maintenant de considerer des /(Mi,M 2 ) bien choisis. 

Remarquons que Ji est lineaire, et que les regies de derivations en chaine s'appliquent si / est 
un produit. II suffit done de calculer Ji pour des fonctions / assez simple et ensuite les combiner. 

Lemme 1.1 Pour toutes matrices A et B, on a: 

k-l 

Ji(AM^B) = J2 tT AM i tr M t l ~ lB (IV.1-7) 

1=0 

k-l 

J^Atr (Mf B)) = ^ tr AM[BM\- 1 - 1 (IV.1-8) 

1=0 

ou, sous forme de series generatrices: 

Jl ( A ^w B ) = tr A ^w tr ~^w B ( IV - 1 " 9 ) 

\ x — Mi J x — Mi x — Mi 

Ji f Atn ( l —B J J = tr A l — B l — (IV.1-10) 

V \x- M 1 J J x - Mi x - Mi y ! 

Les relations eq. ()1V. 1-7)1 , eq. pv. 1-9)1 sont parfois appelees " split-rule" , et les relations eq. (jlV. 1-8)1 , 
eq. ljiV. 1-10)1 sont parfois appelees " merge-rule" . 

Remarque: Les equations de boucles du type eq. (|IV. 1-6)1 . sont a prendre au sens formel, e'est 
a dire qu'elles sont vraies ordre par ordre en puissances de T et en puissances de 1/x lorsque 
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l'on fait intervenir des series geneeratrices. Elles peuvent aussi se demontrer directement par la 
combinatoire des surfaces discretisees [7TJ [72] vue au paragraphe IlEJ Les equations de boucles 
sont un outil standard pour les physiciens de la gravitation quantique EDI 12 UJJ Ulj • Pour le 
modele a deux matrices elles ont d'abbord ete consideree par M. Staudacher jSH], puis systematisees 
par moi meme [Pll], [P9], [P3], [P5], [PI]. 

1.2 Equation de boucles maitresse 

Les equations de boucles donnent des relations entre les moments. II s'agit maintenant de trouver 
des changements de variables /(Mi, M 2 ) bien choisis, afm de trouver une relation fermee entre un 
nombre fini de moments, et ainsi pouvoir calculer tous les moments et observables. 

Definition 1.2 Definissons les observables suivantes (certaines ont deja ete definies au para- 
graphe II \5.Jj\j : 

Wi(x) := I ( tr -\-\ (IV.1-11) 



N \ x - Mi 

Y(x) := V((x) - TW x {x) (IV.1-12) 

^ : =K tr ^> (IVM3) 

X(y) .= Vl{y) ~ TW 2 (y) (IV.1-14) 

U(x , s) : =l/t r _L_ (IV.1-15) 

N \ x — Mi y-M 2 I v ; 

p(,, b) : '( ^M-TO^-W) ) (IV ,. 16) 

TV \ x — Mi y — M 2 I 

R( X ) := L (tr 7 -^Yt{M,)) (IV.1-17) 

[/(x, y; x') := / tr l — YM^MM^l tr —L— \ - N 2 U(x, y)W 1 {x') (IV.1-18) 

Remarquons que U(x,y) et U(x,y; x') sont des polynomes de degre < c?2 — 1 en y, et P(x,y) est 
un polynome de degre cfi — 1 en x et de degre d 2 — 1 en y. On introduit aussi: 

Definition 1.3 Le polynome E(x, y), de degre d% + 1 en x et de degre d 2 + 1 en y, est defini par: 

E(x,y) := (V{(x) - y)(V 2 \y) - x) - TP(x,y) + T (IV.1-19) 

Ce polynome va jouer un role crucial dans toute la suite. La courbe E(x,y) = sera appelee 
"courbe spectrale". 
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Theoreme 1.2 (Eynard [P5], [PI], [P3]j Pour tout x et y, on a (au sens des series formelles): 



(y - Y(x))(TU(x, y) - V 2 \y) + x) + ^U(x, y; x) = E{x, y) 
En particulier, en prenant y = Y(x) on obtient 

Theoreme 1.3 (Eynard [Pll], [P9], [P5], [PI], [P3]) Equation maitresse: 



(IV. 1-20) 



E(x,Y(x)) = ^U(x,Y(x);x) 



(IV. 1-21) 



Preuve: 

L'existence d'une courbe algebrique dans la limite iV grand avait ete prevue par M. Staudacher 
[HE] et explicitee pour des potentiels cubiques. La preuve la plus generale est donnee dans [P3], 
et elle avait deja ete esquissee, dans la limite iV grand, dans [Pll] et [P9]. 

Considerons le changement de variable suivant: 



f(M 1 ,M 2 



V&y) - V±{M 2 ) 



x - M 1 y - M 2 
Le theoreme IV J1.11 avec le Lemme IV J1.1[ donne: 



(IV. 1-22) 



TW 1 (x)U(x,y) + — U(x,y;x) 

= V((x)U(x, y) - P(x, y) - yU(x, y) + V 2 \y)W 1 (x) - R(x) 

et avec le changement de variable suivant pour la matrice M 2 

1 



M 2 M 2 + 5 



x - Mi 



le theoreme IV I1. II donne: 



(IV.1-23) 



(IV. 1-24) 



R(x) = xW x (x) - 1 



(IV.1-25) 



On en deduit: 



T 



(y - V((x) + TWi{x))U(x, y) + j^U(x, y; x) = (V 2 '(y) - x)W 1 (x) - P(x, y) + 1 (IV.1-26) 
Et, en ecrivant Y(x) := V[(x) — TWi(x), on a: 

(y - Y(x))U(x, y) + ^U(x, y; x) = (V 2 \y) - x)W l (x) - P(x, y) + 1 (IV.1-27) 

En multipliant par T, on obtient le theoreme IV 11.2I En corollaire, en prenant y = Y(x), on 
obtient le theoreme IV 11.3I 
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1.3 Fonction de correlation mixte 

De fagon tres similaire, on peut obtenir une equation de boucle pour l'observable mixte 



W h2 (x,y) :=i( 



tr 1 1 



x—Mi y-M 2 

Posons 

WWz,y;j/0 := (tr 1 1 tT -J—- ) - N 2 W 1;2 (x,y)W 2 (y') (IV.1-28) 

\ x — Mi y — M 2 y — M 2 I 

Theoreme 1.4 (Eynard [P5], [PI], [P3]J Wi^O^, y) satis fait V 'equation: 

TW ( \ _ 1 E{x,y) ( U{x,y,x) ^^(x, y; y) 

l ^ V) (x-X(y))(y-Y(x)) + \(x - X{y)){y - Y{x)) x - X{y) 

(IV. 1-29) 

preuve: 

Avec le changement de variable 

M 2 ->M 2 + - + - (IV.1-30) 

2x-M lV -M 2 2y - M 2 x- M 1 v ; 



on obtient: 



TW 2 (y)W 1:2 (x, y) + ^W h2 . >2 (x, y; y) 

= V 2 \y)W h2 {x, y) - U{x, y) - xW lj2 (x, y) + W 2 (y) (IV.1-31) 

i.e. 

(x - V 2 f (y) + TW 2 (y))W^(x, y) + -^W 1>2;2 (x, y; y) = W 2 (y) - U(x, y) (IV.1-32) 

Posons: 

X(y):=V 2 \y)-TW 2 (y) (IV.1-33) 



on a: 



(x-X(y))TW 1 , 2 (x,y)+'^W 1 , 2 . 2 (x,y;y) = V 2 \y) - X(y) - TU(x,y) 

, x E(x, y) — ^U(x, y: x) 
= x-X(y)- { y _Y {x) (i V ' 1 " 34) 



QED. 



Remarque: eq. (jIV.l-29j) dans la limite N grand (i.e. si on neglige le terme en 1/N 2 ), on a: 

qui n'est pas sans rappeler eq. (jIII.12-2|) pour n fini (il faut noter que eq. (|IV.l-29|) a ete calculee 
pour le modele formel, alors que eq. (jIII.12-2|) a ete calculee pour le modele normal). 
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2 Limite N grand 



Notons que pour le modele formel, par definition, les limites N grand et le developpement en 
puissances de 1/N 2 existent. Toutes les fonctions vues precedement se developpent done: 

Wi(x) = Wi°\x) + — wf^x) + ... (IV.2-1) 

: = V((x) - TWi°\x) , Y {1) (x) := -TW[ 1] {x) . . . (IV.2-2) 
E(x,y) = E®(x,y) + ^E^(x,y) + ... (IV.2-3) 

etc... 

Dans la limite N grand, on enleve simplement le membre de droite de eq. (|iV.l-21|) et Ton 
obtient une equation algebrique pour comme fonction de x. On peut egalement calculer la 
resolvante W 2 {0 \y) et X<®(y) = Vtfy) - TW^\y) par la meme methode. On a done: 

Theoreme 2.1 Y^\x) etX<®(y) satis font la meme equation algebrique: 

E {0 \x,Y {0 \x)) = , E w (X {0 \y),y) = (IV.2-4) 

Comme dit plus haut, une relation algebrique etait connue depuis les travaux de M. Staudacher 
[HHl, et l'equation eq. (|lV.2-4|) a ete trouvee explicitement dans [P9] et [Pll]. 

Remarque 2.1 On doit avoir: 

y(0) oX (0) = Id (IV.2-5) 
conformement a ce qui a ete trouve par Matytsin [57] puis prouve par jl^j (voir aussi [ZED- 

Remarque 2.2 Rappelons que W^(x) a ete defmie comme une serie formelle en puissances de 1/x. 
Nous voyons ici que cette serie a un rayon de convergence non nul. C'est une fonction analytique, et plus 
precisement algebrique de x. Elle possede des coupures (et done un rayon de convergence fini). 

2.1 Courbe algebrique 

Nous avons done une courbe algebrique S, d'equation E^ \x,y) = 0, c'est a dire une surface de 
Riemann compacte de genre g, sur laquelle sont definies deux fonctions x et y, telles que pour 
tout point p de la courbe £ on ait: 

Y m (x(p)) =y(p) (IV.2-6) 

La fonction Y^°\x) est solution d'une equation de degre d,2 + 1, elle est multivaluee, elle possede 
d 2 + 1 branches. Autrement dit, la fonction x(p) n'est pas injective, pour chaque x, il existe 
exactement 0I2 + 1 points p 6 £ tels que x{p) = x. 

Commengons par etudier ces branches pour x grand. 
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Points a l'infini 



Theoreme 2.2 Les fonctions x et y possedent seulement deux poles sur £, de diviseurs [x] = 
oo x + d 2 oo y , et [y] = oo y + d 1 oo x , et Von a au voisinage de ces poles: 

Vfc = 1, . . . , c?i , Res^x^ydx = -g k (IV.2-7) 

VA; = 1, . . . , d 2 , ReB O0v y- h xdy = -g k (IV.2-8) 
ReSoo x ydx = Res^xdy = T (IV. 2-9) 

preuve: La resolvante W\{x) est definie comme une serie formelle en puissances de 1/x telle que 
Wx{x) = J + 0{^). II existe done une solution de ea. (llV.2-4D . telle que Y^\x) ~ V({x) - J + 
0(x -2 ), ce qui implique qu'il existe au moins un point oo x G £ tel que: 

x(oo :E ) = oo, y(p) ~ V[{x(p)) - -^T + 0(x(p))- 2 (IV.2-10) 

De meme, pour y grand, il est clair qu'il existe une solution de eq. (|IV.2-4|) . telle que x ~ 

V 2 '(Y^(x))) — Y(o)( x ) + 0(y(°)(x) -2 ), il existe done au moins un point oo^ e £ tel que: 

y(oo tf ) = oc, x{p) ~ V r 2 '( 2 /(p))-^ T + 0( 2 /(p))- 2 (IV.2-11) 

Autrement dit, pres de oo y , y(p) ~ x (p) 1/d2 , et comme il y a d 2 racines c^-iemes de l'unite, il y 
a g?2 branches de la fonction Y(x) pres de oo r Comme la fonction Y^(x) possede exactement 
d 2 + 1 branches, cela implique que les fonctions x(p) et y(p) n'ont pas d'autres poles que oo^ et 
oo^, et que: oo x . est un pole simple de la fonction x(p), et un pole de degre d\ de la fonction y(p), 
et oo^ est un pole simple de la fonction y(p), et un pole de degre d 2 de la fonction x{p). QED. 

Fractions de remplissage 

La definition du modele de matrice formel II I4.5I implique. que pour T assez petit: 

Proposition 2.1 II existe d\d 2 contours disjoints Ai, % = 1, . . . , d\d 2 dans le plan complexe de x, 
tels que: 

Mi = 1, . . . , d x d 2 - 1 & Y(x)dx = 2mTti := 2mr]i (IV.2-12) 

ou les €i sont les fractions de remplissages donnees du modele formel (voir section et oil les 

7]i sont definis par rji := Tei. 

idee de la preuve: en effet, pour T tendant vers 0, les At peuvent etre choisis comme des cercles 
de rayon 0(|T|) et centres sur les d\d 2 valeurs de x telles que V 2 (V{(x)) — x = (voir section 
110]). L'integrale de la resolvante W\(x)dx sur un tel contour est le "nombre de valeurs propres" 
de Mi a l'interieur du domaine delimite par ce contour, e'est done par definition —2me{. Dans la 
limite N — > oo, la resolvante et la courbe £ sont des fonctions analytiques de T (avec un rayon de 
convergence T c ), et Ton obtient des contours Ai pour T fini par continuity. 
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Cycles non-triviaux 

Les contours Ai peuvent etre releves en cycles non-triviaux, notes aussi Ai sur la courbe algebrique 
£. 

Sur une surface de genre g, il y a 2g cylces non-triviaux independants, et il est possible d'en 
choisir une base canonique At, Bi, i = 1, . . . , g, tels que: 

Ai n Aj = o , Bi n Bj = o , Ai n B, = <J<j (iv.2-13) 

Les contours satisfont: 

Vz = 1, . . . , did 2 - 1 f Y(x)dx = 2mTei := 2mrn (IV.2-14) 

Genre 

Une courbe algebrique de la forme eq. (|lV.2-4|) avec E de la forme eq. (|lV.l-19|) a un genre g < 
d\d2 — 1, et a generiquement le genre maximal. 

Une courbe de genre g possede exactement g cycles disjoints non-triviaux independents. 
Autrement dit, la courbe ne peut etre degeneree (i.e. g < d±d 2 — 1) que si certaines combinaisons 
lineaires entieres d'integrales eq. (|IV.2-14J) sont nulles, autrement dit si certaines combinaisons 
lineaires entieres des fractions de remplissage sont nulles (pour T petit, si certaines fractions de 
remplissages sont nulles). 

Nous supposerons a partir d'ici, que E^°'(x, y) est un polynome connu, de la forme eq. (jiV.l-19j) . 
et que les conditions eq. fjl V.2-9jl . eq. fjl V.2-7J1 . eq. fjl V.2-8jl et eq. fjl V.2-14jl sont satisfaites. 

Feuillets 

La fonction x{p) n'est pas injective, pour chaque x, il existe d<i + 1 points p tels que x{p) = x. On 
note: 

Definition 2.1 

x(p) = x <-> p e (p (x),p 1 (x), . . . ,Pda(x)) (IV.2-15) 

et Von suppose que Po(x) est tel que po(x) — * oo x lorsque x — ► oo ; et pour k > 1 Pk{x) — > oo^. 
On decoupe ainsi la courbe £ en d 2 + 1 domaines disjoints appeles "feuillets" , £ , . . . ,£d 2 , tels que 
Pk{x) G £k, et la fonction x{p) restreinte a £k est une bijection. 

Le feuillet numero 0, i.e. celui qui contient oo x est appele feuillet physique. C'est seulement 
dans ce feuillet que la fonction Wi(x) = if(V{(x) — y(po(x))) se developpe en puissances de 1/x 
pour x grand, et coincide avec la fonction generatrice des moments eq. (jll.5-ll|) . 

De meme on note: 

y(p) = y «-> p e (p (y),pi(y),...,p dl (y)) , Ph(y) e £ k (IV.2-16) 
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Le feuillet numero 0, i.e. celui qui contient oo y est appele feuillet physique. C'est dans ce feuillet 
que la fonction W 2 {y) = hiyiiv) ~ x (Po(y))) se developpe en puissances de 1/y et coincide avec 
la fonction generatrice des moments eq. (jll.5-12|) . 

Le choix des domaines £j~ (resp. et de leurs frontieres est essentiellement arbitraire. Les 
images des frontieres par x(p) (resp. y{p)) sont appelees les coupures. 

La seule contrainte est que les points multiples de x{p) (resp. y{p)) soient sur les coupures, ce 
sont les points de branchements. 

Points de branchements et coupures 

Les points multiples de x(p) sont tels que 3k ^ I tels que pk(x) = pi(x), ce sont les zeros de la 
differentielle dx{p). La differentielle dx(p) possede deux poles: un pole double en oo x et un pole de 
degre d 2 + 1 en oo r cela implique qu'elle doit avoir d 2 + 1 + 2g zeros, que Ton note e\, . . . , ed 2 +i+ 2g : 

dx(p) =0 <-> p e (ei, . . . , e d2+1+2g ) (IV.2-17) 

De meme la differentielle dy{p) possede deux poles: un pole double en oo y et un pole de degre 
d\ + 1 en oo x . cela implique qu'elle doit avoir di + 1 + 2g zeros, que Ton note: 

dy{p) =0 ^ pe (ei, . . . , e dl+l+2g ) (IV.2-18) 

Definition 2.2 Les zeros ex, ... , ed 2 +i+2g (resp. e%, . . . , e^+^g) de la differentielle dx (resp. dy) 
sont appeles "points de branchements" de x (resp. y). 

Supposons ici que tous ces zeros sont simples et distincts. Pres d'un zero simple e de dx, la 
fonction Y(x) se comporte en racine carree; 

Y(x) - Y(e) ~ ^x - x(e) (IV.2-19) 

Remarque 2.3 Points critiques 

Les cas ou les zeros de dx et dy ne sont pas simples et distincts, s'appellent points critiques. L'etude 
des points critiques est cruciale dans le cadre des applications des modeles de matrices aux theories 
conformes, mais sort largement du cadre de cette habilitation. Voir par exemple |521 1161 12T1 11UI I12j ... 

Notons que si Ton est a un point critique, et que e est a la fois un zero de degre q — 1 de dx et un 
zero de degre p — 1 de dy, on a une singularity rationelle: 

Y(x) - Y(e) ~ (x - x{e)) p/q (IV.2-20) 

Les points critiques decrivent done des surfaces de Riemann singulieres. II est connu que les modeles de 
matrices avec singularites rationelles d'exposant p/q, decrivent les modeles minimaux conformes (p,q), 
voir [THj . 

2.2 Elements de geometrie algebrique 

Voir jSHJ EZ1 EES] pour une introduction. 
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differentielles holomorphes 

Sur une surface de genre g, il existe g differentielles holomorphes (i.e. sans poles) independantes 
dui(p), i = 1, . . . , g. On peut les choisir telles que: 

dmip) = 5 k>l (IV.2-21) 

On definit alors la matrice des periodes par les integrates sur les £>-cyles: 

r k>l := I d Ul (p) (IV.2-22) 

Un resultat classique de geometrie algebrique [SZ1EH1 montre que r k ,i = et Imr^ > 0. 
Fonction d'Abel 

Donnons nous un point po G S, arbitraire. On definit alors la fonction d'Abel: 

SMUiAiUiBi) -> C a 

rp 

p — > u(p) oh Ui(p) := / dui (IV. 2-23) 

ou le chemin d'integration ne doit croiser aucun des cycles Ai, Bj. 

Cette fonction s'etend naturellement comme un plongement de S dans la Jacobienne C 9 / (Z 9 + 
tZ 9 ): 

8 -> C 9 /(Z 9 + tZ 9 ) 

p -> [u(p)] modulo Z 9 + tZ 9 (IV.2-24) 

Fonctions theta 

Etant donne une matrice r de taille g x g, telle que Imr > 0, on definit la fonction entiere suivante: 

Qg > Q 

9: u — ► 6( U ,r) = j: neZg e^ nt ™)e 2 ^ nt ^ (IV.2-25) 

Par abus de notation, on ometra souvent la dependance en r. 

La fonction 9 possede les proprietes suivantes, pour tout n G Z 9 : 

d(u + n) = 9(u) , d(u + rn) = 9(u) e ~^^)+^)) , Q(-u) = 9(u) (IV.2-26) 

L'ensemble des zeros de la fonction 9 est une sous-variete de codimension 1 de C 9 , notee [9]. 
En particulier, il est clair que les points demi-entiers suivants appartiennent a [9]: 

n + rrn 



2 

lis sont appeles demi-periodes impaires. 



e [9] si (m*n) est impair (IV.2-27) 
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Formes premieres et fonctions Theta 

Sort z une demi-periode impaire. Considerons la forme differentielle holomorphe suivante: 

9 d9(u) 



d Ui (p) (IV.2-28) 



Alors, on definit la forme primaire sur £ x £ 



EM -- $ -^mrwTT (IV - 2 - 29) 

dont on peut montrer qu'elle ne depend pas du choix de z, ni du choix de l'origine po- E(p,q) 
s'annule si et seulement si p = q. 

Noyeau de Bergmann 

Le noyeau de Bergmann B(p, q) est l'unique forme differentielle bilineaire, qui, vue comme forme 
differentielle en p possede un unique pole double hp = q, sans residu, tel que dans n'importequelle 
parametrisation x{p) on ait: 

B(p, q ,~ ^ )d fl + fini (IV.2-30) 
p^q \x(p) - x{q)y 

et telle que: 

Vi = l,...,g <f B(p,q) = (IV.2-31) 

JpeAi 

Cette forme est unique, elle est symetrique B(p,q) = B(q,p), et elle est donnee par: 

B(p, q) := d p d q In 6(u{p) - u{q) + z) = d p d q In E(p, q) (IV.2-32) 

On a la propriete: 

'i = l,...,g J) B(p,q) = 2iirdUi(q) (IV.2-33) 



Vi 



Connexion projective 

Etant donne une fonction f(p) a valeur complexes, non constante, on definit: 

l s * ): ^(S-») (iv - 2 - 34) 

Differentielles de 3e especes 

Etant donnes deux points p\ et P2 distincts sur £, il existe une unique forme differentielle dS PltP2 (p) 
telle que: dS Pl . P2 possede un pole simple a p = p± avec residu +1, un pole simple a p = P2 avec 
residu —1, et pas d'autres poles, et 

i = l,...,9 <f dS PuP2 (p)=0 (IV.2-35) 
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Cette forme est unique, et donnee par: 



E(p,pt) 



E(p,p 2 ) 



pi 



dS pl , P2 (p) = din „, ( = / #(p,g) (IV.2-36) 



<?=P2 



ou cette derniere integrate est calculee le long d'un chemin qui ne croise aucun cycle Ai, Bj. 
On a la propriete: 



Vi = 1, . . . ,g f dS PltP2 (p) = 2m {u^px) - Ui{p 2 )) (IV.2-37) 
2.3 Observables, energie libre, et leurs derivees 

Les resultats presentes dans ce paragraphe sont classiques et ont ete developpes par de nombreux 
auteurs. Voir [P3], [SSI EH EH Hj ■ Us sont presentes pour rendre l'expose complet, et car ils seront 
utilises au chapitre V. 

Theoreme 2.3 Les fonctions a deux points sont donnees par: 

W$(x(p); x{q)) dx(p) dx(q) = B(p, q) - dx(p)dx(q) (IV .2-38) 

(x(p) — x{q)) z 

W&W); v(q)) dy{p) d y ( q ) = b( p , q ) - ^M^L (iv.2-39) 

W${x(p); y{q)) dx{p) dy(q) = -B(p, q) (IV.2-40) 

ou B est le noyau de Bergmann, et ou les Wj^ sont les limites N ^ oo des observables W^-i qui 
ont ete defmies au paragraphe 11 15.41 

Wi.i(x;x') := ( tr l -— tr - l -— ) (IV.2-41) 

' v ' 1 \ x-Mi x' -Mi/ C v 1 

W 2 . 2 (y- y')-= { tr tr > (IV. 2-42) 

,K \ y~M 2 y'-M 2 / c 1 ; 

W!; 2 (x; y) := / tr — i— tr -^—r) (IV.2-43) 
\ x - Mi y - M 2 / c 

Corollaire 2.1 £fa particulier, a points coincidants: 

W$(x(p);x(p)) = \s x (p) , W^(y(p);y(p)) = l -S y (p) (IV.2-44) 
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Theoreme 2.4 Les derivees de la resolvante sont donnees par: 



dTWf\x{p)) 



dgk 



x(p) 



dx(p) = x(p) dx(p) + — Res x(q) B(p,q) 

k q^oox 



dTWf\x(p)) 



dgk 



x(p) 



dx{p) = -- Res y{q) k B(p,q) 



dTW[°\x(p)) 



dT 

dTW[°\x(p)) 



dx(p) = -dSoo^ip) 



x(p) 



drji 



dx(p) = —2indui(p) 



x(p) 



2.4 Energie libre et ses derivees premieres et secondes 

De meme, les formules suivantes sont classiques (voir [P3], [H EH EU Hi] ■ 



Theoreme 2.5 La limite N grand de I 'energie libre est donnee par: 



\fc=0 y/c fc=0 3K i=l 



dF (o) ^ dF (o) dF (o) 1 

— h T— Res xy dx 

ar\i dT 2 oo x 



avec: 



dF(0) 1 T, k J 

— — = — Resx ydx 



dF (o) i 9F (o) 

— — = 7 Resy zrtfy , — — 
og k k oo y drji 



ydx 



(IV.2-45) 

(IV.2-46) 
(IV.2-47) 
(IV.2-48) 



(IV. 2-49) 



(IV.2-50) 



dF<® 
dT 



P (y - V{{x) + -)dx + [ P (x - Vi{y) + -)dy 

X J OOy y 

+Vi(x(p)) + V 2 (y(p)) - x(p)y(p) - T\nx(p) - T\ny(p) (IV.2-51) 



cette derniere expression etant independante du point p. On pent la voir comme une regularisation 
de Vintegrale f°° y ydx. 

Remarque 2.4 si Ton introduit un coefficent h devant le TVM1M2: 

Z - I dM 1 dM 2 e~^ Tr [Vi(Mi)+%(M2)-/»AfiM a ] 



alors on a: 



— - Res xy dx 

2 OO x 



dF (0) 



Oh 



(IV.2-52) 



(IV.2-53) 



h=l 



ce qui rend l'expression ea lTV.2-49l homo gene. 
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Theoreme 2.6 Les de.rive.es secondes de F^ ' sont donnees par. 



dgkdgj 
dgkdgj 



ReB 00a Bes 00w x{p) k x(q) i B(j,, q) (IV.2-54) 
Res 0Ow Res OOs y(p)*y( g yS(p, g) (IV.2-55) 



^2 F (0) 

— — = Res^Res^x^^g)^, g) (IV.2-56) 
ugkogj 

Q 2 _p(0) ^2 _p(0) 

— - = -2m r y , = -2i7r(Wi(oo x ) - Ui(oo v )) (IV.2-57) 

— — = -11177, 7:= bm x(p)— ; r, 7:= hm ^ (IV.2-58) 

...etc. 

2.5 Autres derivees 

Les derivees troisiemes de l'energie libre ont ete trouvees par plusieurs auteurs jHSl Hj , elles se 
reduisent esentiellement a la formulle variationelle de Rauch jHHESj- Pour calculer les derivees 
suivantes, il suffit en principe d'appliquer recursivement la formule variationelle de Rauch, mais 
en pratique les expressions obtenues ainsi sont assez fastidieuses, et leur structure ne saute pas 
aux yeux. 

Une methode diagrammatique a recement ete inventee, par moi meme [SS| pour le modele a 
une matrice, et a ete generalises au modele a deux matrices avec mon etudiant N. Orantin |H5j . 
Elle permet d'ecrire directement n'importe quelle derivee, comme une somme de diagrammes de 
Feynman en arbres. Cette methode est non recursive, et met en evidence une structure tres riche. 
Helas, elle sort trop du cadre de cette habilitation... 



3 Developpement topologique 

On peut developper l'equation eq. (|IV.l-21|) au premier ordre en 1/N 2 : 

Wi(x) = W[°\x) + —w£\x) + ... (IV.3-1) 

Nous venons de voir que (x) est solution d'une equation algebrique. Nous voulons maintenant 
calculer W^\x). Ceci a ete fait dans [P5](pour g = 0) ou [PI] (pour g quelconque). Je vais 
presenter ici, une methode differente de [P5]ou [PI], et inspiree de mon recent travail [SSI EE] 
pour obtenir le theoreme suivant: 
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Theoreme 3.1 (Eynard, Kokotov, Korotkin [P5], [PI],), (Eynard, Orantin \33$). La resolvante 
a I'ordre 1/N 2 est donnee par: 

TWf\ X ( p) ) MP) = \ T Res Res ,,.,^^^ Wrf » (IV.3-2) 



t q^p>^ q (y(q) - y(p'))(x(q) - x{p')) 



Idee de la preuve: 

L'equation maitresse eq. (jIV.l-21|) developpee au premier ordre en 1/N 2 donne: 

TwV>(r\- T2 U i0) (x,Y(xy,x)+TPW(x,Y(x)) 

1Wl {X) ~ E y (x,Y(x)) [lV - 6 ' 6) 

De plus, on a (nous ne detaillerons pas ici le calcul (voir [P5], [PI]), qui permet de trouver 
U^{x,y;x') a partir de eq. fl V.1-2njl et eq. ()]V.2-3 



T 2 U^(x( q )MQ);x(q)) dx{q) = _f^ B(q,p k (q)) 

E y(x(q),y(q)) ^ (y(q) - y(Pk(q))) dx(p k (q)) 

ou les Pk(q) sont les solutions de x(pk(q)) = x(q) dans les autres feuillets (on suppose po(q) = q). 

En manipulant eq. (jIV.3-3|) . on peut montrer (cf [PI]) que w[ l \x(p)) dx(p) est une different ielle 
sur £. L'hypothese eq. (jlV.2-14|) du modele formel, implique que: 



Vi = 1, . . . ,g & Wl \x) dx = (IV.3-5) 



r(l), 

VY 

'A, 

et implique aussi que W^\x) dx ne possede pas d'autres residus sur £ . En particulier W^\x) dx 
ne possede pas de poles aux zeros de E y (x,Y(x)) autres que les points de branchements. Tout 
ceci implique que W^{x) dx est une different ielle sur £ avec des poles (de degre 4) seulement aux 
points de branchements e^, et qui satisfait eq. (jIV.3-5|) . 

La formule de Cauchy permet d'ecrire (pour tout a G £): 

W[ 1] {x{v)) dx(p) = -Res dS qa (p)wi 1] (x(q))dx(q) (IV.3-6) 

ou le residu est calcule comme une integrate le long d'un petit cercle entourant le point p. On peut 
deformer ce contour d'integration, de fagon a transformer le residu en p, en residus aux points 
de branchements. En vertu de l'identite bilineaire de Riemann (HBl EZj, et grace aux relations 
eq. (jIV.2-35j) et eq. (|IV.3-5|) . il n'y a pas de contribution des cycles A%, Bj. II reste done: 

d 2 +l+2g 

W[ l) \x(p)) dx(p) = V Res dS qa (p)wj 1] (x(q))dx(q) (IV.3-7) 

i=l 

On insere done eq. (jIV.3-3|) dans eq. (jIV.3-7j) . et Ton remarque que P E ^y{x)) ^ x n ' a P as ^ e 
residu aux points de branchements car P^(x,y) est un polynome. D'ou: 

W^(x(p)) dx{p) = -T £ 9 Res dS q M dx{q ) (iy.3-8) 
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Et done, en utilisant eq. (jIV.3-4|) : 

d 2 +l+2g d 2 



E Jg^ >0 (p) (y(g) _ y(p /))^(p - ^(g)) 

_ V Res ^9 Res 

tT ™ ^ " Kp'))W) " x(q)) 

d 2 +l+2g di « 

- E S^^E^iFS^i) 

d2+1+2 9 5("g 

= Res dS qjCt (p) Res 



+ Y Res dS a (p) y B ^Mq)) 
^ Res Res B M 



g ^e lP >-+g (y(q) - y(p'))(x(q) - x(p')) 
dSq !P >{p)B(q,p') 



- ^y Res Res 



(IV.3 - 9) 



^ q-*e il y-, q (y(q) - y(p'))(x(q) - x(p')) 



L'egalite entre la ligne 1 et 2 est une simple reecriture de x(pk(q)) = x(q) sous forme de residus. 
L'egalite entre la ligne 2 et 3 est obtenue en deplacant le contour, pour prendre tous les autres 
residus, i.e. p' = q et les d\ autres solutions de y(p') = y(q), que Ton appelle pk(q), k = 1, . . . , d\. 
Dans la quatrieme ligne, on calcule explicitement les residus en p' = pk(q), et Ton voit qu'ils 
n'ont pas de poles lorsque q = Le passage de la cinquieme a la sixieme ligne, est obtenu en 
echangeant l'ordre des residus, il montre que cette expression est bien independante du choix du 
point a. On trouve done eq. (|IV.3-2J) . QED. 

a partir de la, on obtient: 

Theoreme 3.2 (Eynard, Kokotov, Korotkin [P5], [PI] j L'energie libre a l'ordre 1/N 2 est donnee 
par: 

= ~ In (~gXl[ d2 rf ''ffVte)) (IV.3-10) 



24 

ou 



Ci(p) := y/x(p) - x{e t ) (IV.3-11) 
y'iet) := (IV.3-12) 
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et t x est determinee par [PI], [5$ : 



La formule variationnelle de Rauch implique que eq. ()IV.3-13|) definit bien une fonction r x , 
parfois appelee fonction r de Bergmann [53]. On peut definir une fonction r y de fagon similaire, 
et on peut verifier que eq. (|3.2|) est bien symmetrique dans l'echange x <-> y. 
Idee de la preuve: 

II s'agit de verifier que l'expression eq. (jlV.3-10|) satisfait bien toutes les relations du type 
eq. ()H.5-2j) . i.e. par exemple: 

— — = -- Resx fc Wr\x)dx (IV.3-14) 
dg k k oo x 

On calcule le membre de gauche de eq. (jIV.3-14|) . i.e. les derivees par rapport a g^, en utilisant 
les relations du theoreme IV I2.41 et on calcule le membre de droite de de eq. (jlV.3-14j) en utilisant 
le thoreme IV 13. II La preuve, fastidieuse. est exposee dans [PI], et elle est trop longue pour 
aparaitre ici. 

Par ailleurs, la fonction t x est connue pour etre proportioned au determinant d'un Laplacien 
sur 8 (voir [PI], PI EE!)- 

Le developpement systematique de la resolvante et de ses derivees a ete obtenu tres recement 
par (Eynard, Orantin) [SH], par une methode diagramatique, que je ne presenterai pas ici. Toute- 
fois, cette methode diagrammatique n'a pas encore pu s'appliquer au calcul du developpement 
de l'energie libre en 1/N 2 . Ceci reste a comprendre. Pour l'instant, il n'existe pas d'expression 
connue pour par exemple. 
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Chapitre 5 

Asymptotiques des polynomes 
biorthogonaux 

Dans ce chapitre, nous allons voir comment on peut deviner les asymptotiques a n grand, pour 
les polynomes biorthogonaux. Les chapitres precedents n'etaient qu'une introduction a ce calcul, 
pour poser les notations, et mettre en place les elements necessaires. 

L'idee est la suivante: la formule de Heine eq. (jIII.l-2|) permet de voir les polynomes orthogo- 
naux comme des fonctions de partition d'un modele normal. Le modele normal est lui meme la 
transformee de Fourrier d'un modele normal a symmetrie brisee. Le modele a symmetric brisee, 
a, sous certaines hypotheses, le meme developpement asymptotique qu'un modele formel. Par 
la methode des equations de boucles, nous pouvons trouver les asymptotiques du modele formel 
en termes de geometrie algebrique. II suffit alors de refaire la transformee Fourrier inverse pour 
obtenir les asymptotiques des polynomes. Cett methode a ete introduite dans [Pll]et [P9], et a 
ete presentee aux Houches en 2004 [31 j . 

Les asymptotiques ainsi obtenus s'expriment naturellement en termes de geometrie algebrique 
d'une courbe complexe, et Ton voit aisement que cette courbe complexe n'est rien d'autre que la 
limite n grand de la courbe spectrale vue en eq. (jIII.7-9|) . 

II faut noter que les asymptotiques presentes ci-dessous ne sont que des conjectures (a ce 
jour), soutenues par des arguments physiques tres raisonables, mais ils ne sont pas demontres 
mathematiquement. La methode la plus prometeuse pour les demontrer est la methode de 
Riemann-Hilbert introduite par (TUJ HI] et [THJ I2DI EI] • H aparait egalement probable que 
Ton puisse un jour justifier rigoureusement les hypotheses presentees ci-dessous, sur lesquelles se 
base l'intuition physique, sans passer par Riemann-Hilbert. 

Les asymptotiques presentes ci-dessous ont ete introduits pour la premiere fois en 1997 dans 
[Plljet [P9]pour le cas ou la courbe spectrale est de genre 0, et ils ont ete presentes a la conference 
de l'AMS a Montreal en 2002 0. 
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1 La courbe spectrale 



Supposons qu'il existe une courbe algebrique S adequation E(x,y) = avec: 

E(x, y) = (V((x) - y)(Vi(y) - x) - P(x, y) + 1 (V.l-1) 
ou degP = (di — 1, d 2 — 1), et telle que pour tout cycle sur S, on ait: 

Re<pydx = (V.l-2) 



Supposons que la courbe est de genre g < d±d 2 — 1, et que Ton ait pu choisir des cycles Ai, 
entourant des coupures [027-1,02/], tels que les: 

e/ := ^— f ydx , e g+l := 1 - J^e/ (V.1-3) 

soient tous positifs. 



A! 



1=1 



Remarque 1.1 L'existence d'une telle courbe est obtenue a partir de la courbe algebrique du modele 
formel, vue au chapitre precedent. Pour chaque e/ > 0, il existe une courbe algebrique, et la partie reelle 
de Penergie libre associee est une fonction convexe des e (en effet, sa derivee seconde est — 2i7rr, dont la 
partie reelle est toujours positive). Elle admet done un unique minimum, et e'est ce minimum que l'on 
considere ici. Le minimum peut etre atteint aux bords, dans ce cas certains e/ sont nuls, et le genre est 
donne par le nombre de fractions de remplissage non nulles. Au minimum d'un e/ / 0, on a ^ff- = 0, 
i.e. Re § B ydx = 0. 



1.1 Feuillets 

Comme nous l'avons vu au chapitre precedent, cette courbe possede deux points a l'infini, notes 
oo x et oo yi et possede di + 1 (resp. d\ + 1) feuillets en x (resp. en y). C'est a dire que pour tout 
x E C (resp. y e C), il existe d 2 + 1 (resp. d\ + 1) points de £ tels que: 

Vi = 0, . . . , d 2 , x(pi(x)) = x , (resp. \/i = 0, . . . , d 1: y(Pi(y)) = y ) (V.l-4) 

On note po(x) (resp. po{y)) le point qui est dans le feuillet physique, i.e. qui s'approche de 00 x 
(resp. oo y ) lorsque x — > 00 (resp. y — > 00). Les autres points Pi(x) (resp. Pi{y)) avec i > 
s'approchent de oo y (resp. oo x ) lorsque x — > 00 (resp. y — > 00). 

1.2 coupures, points de branchelents et bases d'homologies 

Les points de branchement en x (resp. en y) sont les zeros de dx (resp. dy), on les note e$ (resp. 
e?j). On supposera ici qu'ils sont simples et tous distincts. On notera: 

di := x(ei) , (resp. 6i := yfe) ) (V.l-5) 
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Chaque point de branchement du feuillet physique a, (resp. bi), est tel que 3j > tel que: 



Pj(ai) = po(at) , ( res P- Pj( b i) = Po(h) ) (V.l-6) 
Definissons le potentiel effect if: 

rpo{x) rpo(y) 

Vi, eS (i)(x) := / ydx , (resp. V2,eff(i) (j/) := / xdy) (V.1-7) 

•'Pi (as) Jpjiv) 

On a: 

K,eff(i)(ai) = , (resp. V 2 ^ m {h) = 0) (V.1-8) 

et Vi tC ff(i)(x) (resp. V2, e ff(i)(z/)) se comporte en (x — a^) 3 / 2 (resp. (y — 6j) 3 ^ 2 ) au voisinage de a{ 
(resp. bi). 11 existe des directions ou Vi^^x) (resp. V^eff^d/)) est reel el strictement croissant. 
Ceci definit (au moins) un chemin allant de a\ (resp. b{) a 1'infini. Pour chaque coupure /, 
i.e. reliant deuxpoints de branchements [a 2 /-i,ci2/] (resp. [621-1, ^2/]), on a done deux chemins, 
partant chacun d'une extremite de la coupure. Notons T x i (resp. T v i) le chemin allant de 00 
a 00, constitue de la reunion de la coupure et des deux chemins a potentiel effect if s croissant a 
partir des deux bords. Les chemins T x 1 ainsi obtenus sont necessairement sans intersections, et 
vont necessairement a 1'infini dans des directions ou ReVi(x) > (resp. ReV^y) > 0). On peut 
aussi facilement montrer qu'ils sont homologiquement non nuls, et que les T x 1 X T y i forment bien 
une base d'homologie (cf chemins similaires dans p'6\ ) . 

Nous utiliserons cette base pour definir le produit scalaire des polynomes biorthogonaux, i.e. 
il faut multiplier la matrice k vue en III.2-4I qui a servi a definir le produit scalaire des polynomes 
biorthogonaux. Nous supposerons desromais que k est la matrice des coefficients qui decrit le 
chemin T des polynomes biorthogonaux dans cette base: 

r = 53«jfjxfj (v.1-9) 
1 

ir n (x)a m (y)e- (Vl{x)+V2{y) - xy) dxdy = hj nm (V.1-10) 



r 



Nous noterons: 



Kj = e 2twui (V.l-11) 



Remarquons que dans les applications les plus frequentes, les Ki valent ±1, et done les v% sont 
des entiers ou demi-entiers. 

Si un coefficient kj est nul, il faut choisir ej = 0. Le genre de la courbe est donne par le nombre 
de coefficients non nuls. Remarquons que si le chemin est l'axe reel, celui ci est somme d'au plus 
la moitie des chemins, et done le genre est au plus: 

di + 1 do + 1 

9 < (V.l-12) 
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1.3 Notations 

Rappelons que 

Resydx = Res xdy =1 , VJ , — <£ ydx — ej > , — d> ydx = (j (g R) 



VX , Res —^—ydx = V{(X) , VF , Res —^—rxdy = V£{Y) 

X — X oo y y — Y 



(V.l-13) 

On choisit une base canonique de differentielles holomorphes: 

duj := 5 U (V.l-14) 



et la matrice des periodes associee r (on a r* = r et Imr > 0): 



/ duj:=Tu (V.l-15) 



On choisit une demi periode impaire z = n+ ™ avec Ylii n i m i e 2Z+ 1, et on considere 6^(w, r) : = 
#(w + z, r). On definit la differentielle holomorphe: 



(V.l-16) 



v=0 



On definit (sur un domaine fondamental du recouvrement universel, et avec une coupure le 
long d'une ligne [oo x , oo y ]): 

T{p) := V 1 (x(p)) - In (x(p)) + [" (y(q) - V[{x{q)) + dx(q) (V.l-17) 



q=oo x 



T(p) := V 2 (y(p)) - In (y(p)) 4 / ' ( x(q) - V^y{q)) + ) Mq) (V.l-18) 



x(q) 
1 



'q=oo y \ 

Notons que e T ^ n'a pas de coupure le long de [oo x , oo^]. 
On definit: 

H := T{p) + f{p) - x(p)y(p) (V.l-19) 

qui est independant du point p (en effet d\i = ydx + xdy — dixy) = 0). 
On definit aussi: 

./ X 0z(u(p) — ll(00y), r) X(p) _ ,. /NWS ™\ 

Hp) ■= Trrri — > ^ := hm ttt ' ^ := hm a (p v.1-20 

Vziuip) - Ui00 x ), T) P^oo x A{p) p-^ooy 



ff( n ^(P) 0z(u(ooy) -u(oo x )) ~ _ „ dh z jp) 9 z iujoo x ) - ujooy)) 

{Ph 7 cfa(p) 0,(u(p) - «(oo x ))* ' W ' 7 dy(p) 9 z iuip) - u(oo y ))* 1 •" J 
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Notons que H(oo x ) = 1, resp. H(oo y ) = 1, et 

H{p) 7 KF) dx(p) K } 

Posons: 

Vk := N{re -Q-v + k{u(oo x ) - u(oo y )) (V.l-23) 

et 

^ fe - = 2n / 27T77 c _ JV/J e _2 i7 riV6( u (cx) ;c )-«(oo y )) 0(Vk-i,r) , y 1 _ 24 , 

Definition 1.1 Aons definissons les fonction suivantes: 



h{p) ■= X I ^ (7A( P ))- fc e M ^)-^-)) %fc + "ff ~ " (00 - } ' T) (V.l-25) 



/ fc (p) ■= J^M ( JL^| _ c -2i« N eM P )-u(oo v ) AVk-u(p)+u(0O y ),T) 

9u(p) = ±- ^ZJm ( e*^))!!^^^ (V.l-27) 

h(p) ■= ^ A-^(P) (7A(p)r fc e-^W-^-)) gfa - fc - "(P) +"(°°»)^) (v.1-28) 

Notons que fk{p)e~ NT ^ n'a pas de discontinuity le long des cycles £>/, et est multipliee par 
une phase e 2t7T ( u i+ z i) = e ^^i Kj lorsqu'on traverse un cycle Aj. 
fkip) se comporte asymptotiquement a l'infini comme: 



fk(p)e- NT{p) ~ x(p) N - k e- NVl( - x V»(l + 0(l/x)) (V.l-29) 

Definition 1.2 Definissons les matrices carrees suivantes de taille d2 + 1 

T(x) := dmg(T(p (x)),T( Pl (x)),...,T(p d2 (x))) (V.l-30) 

^•(^) := /i(Pi(a:)) , i,j = 0,...,d 2 (V.l-31) 

Gij(x) := <7i(pj(a;)) , i,j = 0,...,d 2 (V.l-32) 
ei Zes matrices carrees suivantes de taille d± + 1 

T(y) := diag(f(p (!/)),T , (p 1 (y)),...,T'(p dl (y))) (V.l-33) 

- /ife(y)) , J = 0, • • • , di (V.l-34) 

:= 9i(Pj(y)) , i, j = 0, . . . , di (V.l-35) 
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2 Conjecture 



Conjecture 2.1 (Eynard [Pll], [P9]j ^ es asymptotiques des systemes fondamentaux des 
polynomes biorthogonaux ( definis par le theoreme \8.1\) . dans le regime 



sont donnes par: 



n^oo,N^oo, \n-N\ =0(1) (V.2-1) 

9 N (x) ~ F(x) e~ NT ^ C x (l + 0(1 /N)) (V.2-2) 

* N (y) ~ F(y) e~ N *M Cy (l + 0(1/ N)) (V.2-3) 

*n(v) ~ G?(y) C y (l + 0(l/iV)) (V.2-4) 

$ N (x) ~ 4(1 + O(l/A0) (V.2-5) 

Cz, Cy, Cj, soni rfes matrices carrees constantes par morceaux. 

Remarque 2.1 Ces asymptotiques ont ete introduits pour la premiere fois en 1997 dans [Pll] et [P9] 
pour le cas ou la courbe spectrale est de genre 0, et ils ont ete presentes a la conference de l'AMS a 
Montreal en 2002 Nj]. Ces asymptotiques sont des ansatz, qu'il faudrait maintenant prouver, par exemple 
par la methode de Riemann-Hilbert ^U] et |18j . 

Remarque 2.2 Les polynomes orthogonaux ^Ar_j(x) sont obtenus dans la premiere colonne de la 
matrice ^^(x): 

i/) N -i(x) = (*N(x)) m (V.2-6) 

Remarque 2.3 Ces asymptotiques sont uniforme pour tout x, hors du voisinage des points de branche- 
ments. Les lignes de discontinuite des C x , sont donnees par: 

rpj{x) 

Re / ydx = (V.2-7) 

Jpi(x) 

Remarque 2.4 La methode euristique pour obtenir ces asymptotiques est presentee dans [Pll] et 
[P9], de meme que dans |SU pour le genre zero. Lorsque la courbe a un genre plus eleve, la sommation 
sur les fractions de remplissages (cf eq lll.3-10|) . produit une fonction 8, comme dans ^3]- De plus, 
cette methode euristique fonctionne a x fixe, i.e. on peut calculer la matrice C x tres explicitement pour 
un x donne, sans avoir a se preocupper des discontinuites. La matrice C x provient de la matrice de 
changement de base d'homologie pour calculer les transformees de Fourrier et Hilbert, par la methode du 
col (methode du col habituelle en dimension 1). Les discontinuites sont obtenues a posteriori. Les lignes 
de discontinuites apparaissent d'abord car les Pj(x) sont par definition discontinus le long des coupures, 
puis, par un phenomene de Stokes lorsque Ton calcule les transformees de Fourrier et transformees de 
Hilbert par la methode du col. Les matrices de Stokes et de saut le long des coupures pourraient etre 
calculees tres explicitement par cette methode euristique. 
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3 Verifications 



Verifions que cet ansatz satisfait certaines proprieties attendues: 



3.1 Asymptotiques 

On a: 



f k {p) ^x( P y k (l + 0(l/x)) (V.3-1) 



fk(p) ~ o(y(p)*"^) (V.3-2) 



ce qui garantit bien que: 



PN-k{x) ~ C.onX^a + Oa/rr)) (V.3-3) 



II faut done choisir C x00 = 1 dans les secteurs ou x — >• oo. 

De fagon generale, on trouve que les systemes fondamentaux donnes dans la conjecture V i2.lt 
satisfont les memes asymptotiques a x grand que ceux donnes au theoreme III I9.11 

3.2 Transformees de Fourrier et Hilbert 

On peut verifier que pour j > 0, 

W-M = [ ^N-i{x) e N *v dx (1 + 0(1/N)) (V.3-4) 



l'integration se fait par la methode du col habituelle. Le chemin d'integration non precise ici, 
depend de x de fagon constante par morceaux, car il depend de la matrice C x , et de la classe 
d'homologie des chemins passant par les points cols. 

On peut aussi verifier par la methode du col et la formule de Cauchy, que 

$nUv)= f f -^ipN-i(x)e Nxv ' dxdy' (1 + 0(1/N)) (V.3-5) 

J J y-y 

3.3 Orthogonalite 

On peut verifier que: 

' ipN-iix) fa-jiv) e Nxy dx dy = 5 tj + 0(1/N) (V.3-6) 




en effet, l'integration se reduit a un residu au point oo^. (ou au point oo y ) sur la courbe algebrique, 
de la forme: 

Jk Reg dh z (p) 9 z (u(oo x ) - ujooy)) j^Ajp))^ 1 Ojrji + ujp) - ujoo x )) flfa - ujp) + ujoo y )) _ 
h P ^oo x 9 z {u{p) - u(oa x ))6 z (u(p) - u(oo y )) 9{rj i )9{rj j ^ 1 ) 13 
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(V.3-7) 

si i > j, il n'y a pas de pole en oo x et le residue s'annule, de meme, si j > i, il n'y a pas de pole 
en oo y et le residue s'annule, et si i — j, le residue vaut bien 1. Ceci montre que le coefficient h n 
(defillED vaut: 

h n ~h N - n (l + 0(l/N)) (V.3-8) 
de meme, le coefficient 7„ (def JIII.2^8]) vaut: 



7 " ~ ^ 77 y (1 + 0(1/JV)) (V.3-9) 

3.4 Systemes different iels, courbe spectrale et dualite spectrale 

Les 4 systemes de la conjecture V I2.1l satisfont 4 systemes d'equations differentiels: 

V x {x) ~ F{x) T'(x) F(x)- 1 + 0(1/N) (V.3-10) 

~ G(x) T'{x) Gix)- 1 + 0(1/N) (V.3-11) 

P 2 (y) ~ f'(y) F( y y l + 0(1/JV) (V.3-12) 

V 2 (y) ~ f (y) G(y)- 1 + 0(1/JV) (V.3-13) 

Les valeurs propres de sont done celles de T'(x), i.e. ce sont les y(pi{x)), autrement dit la 

courbe algebrique que nous avons introduite est bien la limite de la courbe spectrale, et ce pour 
les 4 systemes. Nous retrouvons le thoreme de dualite spectrale: 

Theoreme 3.1 Les 4 systemes differentiels ci-dessus ont la meme courbe spectrale, qui est la 
courbe algebrique £. 

Remarque 3.1 La matrice T>\{x) n'a pas necessairement de limite lorsque ./V — > oo, car F{x) depend 
explicitement de N. 

3.5 Dualite et Christoffel— Darboux 

Theoreme 3.2 La matrice 

A = (F(x)G(x)') -1 (V.3-14) 

ne depend pas de x, et a la forme de la matrice de Christoffel Darboux (cf theoreme III \4 ■ 1[ ecrite 
dans une autre base): Ay — si (i — et j ^ 0) on si (j — et i ^ 0) on si (i > 0, j > 0, et 
i+j> d 2 ). On a A 00 = -Jn-i- 

idee de la preuve: Chacun des A$j peut etre ecrit comme des combinaisons de fonctions sur 
la courbe £ . On verifie que cette fonction est bien une fonction meromorphe sur 8 , i.e. qu'elle 
reprend bien sa valeur apres un tour autour de chaque Ai et chaque Bj, et qu'elle n'a pas de poles. 
Toute fonction meromorphe sans pole etant constante, on trouve que A ne depend pas de x. La 
valeur des A^ peut done etre calculee pour x grand, et Ton trouve le resultat annonce. 
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4 Idee euristique sous-tendant la conjecture 



Nous allons exposer rapidement les idees qui conduisent a cette conjecture. Elles se basent sur les 
articles [Plljet [P9], et sur [T3] (voir aussi [21]). 

4.1 Formule de Heine et modele normal avec potentiel rationel 

La formule de Heine eq. (jIII.l-2|) permet de voir les polynomes orthogonaux comme des fonctions 
de partition d'un modele normal: 

Z NoTm (N, T, r, K) 
Knit) = ^ (at rp n \ (V.4-1) 

Z Norm (7V, 1,0, k) 

ou T = § et r = et: 

Z Kana (N,T,r,K) := / dM 1 dM 2 e~^ * [WH^-mo+v^-m^] = e -^»»Wv) 

(V.4-2) 

autrement dit, on a simplement ajoute un terme logarithmique au potentiel V\, et une temperature 
T = jj. Jusqu'ici, nous n'avons considere que des potentiels polynomiaux, mais il est connu jH] que 
tout les raisonements precedents se generalisent bien au cas ou V{ est une fraction rationelle. De 
plus, ici, le terme logarithmique est multiplie par r = —1/N qui est petit, et nous allons effectuer 
un developpement de Taylor au voisinage d'un potentiel polynomial. 

4.2 Transformee de Fourrier d'un modele normal brise 

Pour un chemin d'integration T = KjY x i x T y i, on a, d'apres eq. (jll.3-10|) : 

Z Norm (N,T,r,n) = ^ nf 1 Z Normh (N, T, r, mi) (V.4-3) 

E/ mi=n 



avec: 



^Normb(^V, T,r,mi) := e" 



-n 2 overT 2 F Normh (N,T,r,mj) 



Y[ / dM xJ dM y 

j J H mi (Y x I )xH mi (Vyi) 



l y,i 

'H mj (p/)xff m; (r9/) 

e -f tr [Vi(M x j)+rln{t-M x j)+V2(M yi r)-M X}I My tI ] 

Y[ det (M xJ ® l ftfcil - l nj ® M M)I ) 
(/)<(fc,0 

JJ det (M yJ <g> l nM - l ni (8) M yjifc J (V.4-4) 
(i)<(M) 

On fait l'hypothese, que si les chemins Ti sont "bien choisis" (i.e. section H~2*|) . alors l'energie 
libre du modele normal brise possede un developpement en puissances de 1/N 2 pour N grand, 
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i.e. on identifie le modele normal brise a un modele formel (a des termes exponent iellement petits 
pres) : 

^Normb(iV, T, r, mj) = Z FoTm (N, T, r, l mj ) = e~£ 'W^*"*) (V 4 _ 5) 

4.3 Developpement a N grand du modele formel 

Nous savons que le modele formel peut se developper en puissances de 1/N 2 : 

- ^ In Z Form (iV, T, r, = F(°)(T, r, m ) + ^F^T, r, ifr) + O(^) (V.4-6) 

ou F^(T,r,r)i), F^ 1 ' (T, r, rji) , etc, sont donnees par Th l2.5l Th l3.21 et sont analytiques dans 
tous leurs parametres. 

Nous poserons F = F^ car nous n'aurons pas besoin de l'ordre suivant. 

En utilisant les theoremes IV 12.41 IV l2.5l et IV 12.61 on a (en ecrivant £ = x(q)): 

dF f 

F v := — = -(b ydx = -2m(i V.4-7 

F r := d -^ = -{T{q)-V l {x{q))) (V.4-8) 
dF 

F T := — =fi (V.4-9) 

F m : = ^g- = -2^ (V.4-10) 



dr- 

d 2 F 



F„:= — = -\nH(q) (V.4-11) 



<9T2 



= -11177 (V.4-12) 



<9 2 F 

F ^ ■= IT^T = -2i7r(«(g) - Moo,)) (V.4-13) 
<9 2 F 

-2z7r(u(oo a; ) -u(oo v )) (V.4-14) 



drjjdT 



Frt:= ^§r =lniMq) (v - 4 " 15) 

Au voisinage de T = 1, r = et 77 = e, en posant m = Nr], la formule de Taylor a l'ordre 2 
donne: 

e -^F(T,r, V ) ^ e -N 2 F(l,0,e) e -N{m-Ne)F v Q NF r e -N(n-N)F T 

e -i {m - Nt) t Fvvim - Ne) e -iF rr e -l(n-N)^F TT 
e Fr, r (m-Ne) Q -{n-N)F vT {m-Ne) Q (n-N)F rT 
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~ ( 7 7)K"- 7V ) 2 e ~N 2 F(l,0,e) e -N{n-N)n 

^/W{^A(qr- N e- N ^- v ^^ 

»7r(m-.!Ve)*T(m--iVe) 2iirJVC*(m-.iVe) 
^-2i-K(u(q)-u(oo x )) t (m-Ne) e 2iir(n-N)(u(oo x )-u(oc y )) t (m-Ne) 

(VA - 16) 

ou toutes les derivees sont calculees en T — 1, r — et r] — e. 
4.4 Sommation sur les fractions de remplissage 

Comme la matrice des periodes tjj a toujours une partie imaginaire positive (voir paragrapheEZ 
ReF(l, 0, if) est une fonction convexe des m. Notons e son minimum. La condition de minimalite 
s'ecrit: 

dBeF 

= = 27rImC/ (V.4-17) 

drji 

ce qui implique que les (i sont reels. La convexite de ReF, implique que dans la somme sur les 
mi, seuls les mi tels que \mi — Nei\ ~ 0(1) contribuent notablement a la somme. ceci entraine, 
d'une part, on peut se contenter de l'approximation de Taylor de l'energie libre au 2e ordre, et 
d'autre part, on peut sommer sure les mi allant de — oo a +oo, ce qui produit une fonction 9: 

Z Norm (iV, T, r, K ) ~ (-yrfW e -N*F(l,0,e) e -N(n-N), 

^Wi^)A(qr- N e- N ^- v ^^ 



v t m^i-K{m-Nt) t T{ra-Nt) 2iirN(? (m— Ne) 



2in(u(q)-u(oo x )Y(m-Ne) 2i7r(n-Af)(u(oo a: )-n(oOj / )) t (m-Afe) 



mi 

e -i... v _ w v ..„ e 

(jj)h( n -W e -N 2 F(l,0,e) e -N(n-N)v 



e i-KNt t rNt e -2inNi; t N 

0(VN-n + u(q) - u(oo x ), r) 



MriVeViVe e -2iwN( t Ne 2m(u(,q)-u(,oo x ))*Ne e -2i7r(n-V)(u(oo I )-u(oo y ))*Afe 



(VA - 18) 

et le rapport Z NoTm (N, T, r, ft)/Z Norm (iV, T, 0, k) vaut dans cette approximation: 
*•„(*(?)) = ^ orm (iV,T,r,K) ^ e -iV(T( g) - 

^NormlA', -T , 0, K) 

e 2wrJVe«(«(«)-«(oo»)) 0(VN-n + u(q) -u(oO x ),t) 

d(VN-n,r) 

(VA - 19) 

i.e. on retrouve la fonction f N _ n (p) definie a l'equation IV. 1-251 

Les arguments presentes ici sont euristiques et indicatifs. lis servent juste a deviner la bonne 
forme a donner a la conjecture. 
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Je pense qu'il doit etre possible de justifier rigoureusement toute cette approche, et done de 
prouver la conjecture, sans passer par une methode de Riemann-Hilbert. Mais ceci reste a faire... 
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Chapitre 6 
Conclusion 



Ce memoire fait un resume de plus de 10 annees de mes travaux sur le modele a deux matrices. 
Les travaux les plus anciens ne concernaient que les courbe spectrale de genre 0, et c'est petit a 
petit que la geometrie algebrique un peu plus evoluee, a fait son apparition, et a permis d'etudier 
les courbes de genre quelconque. 

Mes travaux avaient debute en temps que physicien, le but etant de comprendre un peu mieux la 
gravitation quantique et les theories conformes. Dans ce cadre, les physiciens etaient tres interesses 
par le developpement topologique du modele formel, et seul le cas de genre zero correspondait 
vraiement a la "gravite quantique" (en effet, la fonction de partition est une somme sur des surfaces 
discretisees colees seulement par leur bords, pas par leurs centres). 

L'idee de regarder les courbes de genre plus eleve, c'est imposee naturellement, d'une part par 
souci de completude, d'autre part pour resoudre le "puzzle de Brezin et Deo" (depuis resolu par 
Deift & co PHI); et enfin car la theorie de Dijkgraaf et Vafa a mis en evidence le role joue par les 
courbes de genre plus eleve en theorie des cordes. 

Dans ce memoire, j'ai commence par definir les modeles de matrices, de la fagon dont je les 
comprends, puis j'ai explique ce que Ton sait faire avec la methode des polynomes biorthogonaux, 
et ensuite par la methode des equations de boucles. Le dernier chapitre montre que lorsque Ton 
combine les deux methodes, on obtient une comprehension beaucoup plus grande, et ceci permet 
par exemple de deviner la forme des asymptotiques des polynomes biorthogonaux. Bref, on gagne 
toujours a attaquer un probleme par plusieurs cotes, et a connitre aussi bien les travaux des 
mathematiciens que des physiciens. 

Les perspectives de ces travaux sont nombreuses, tant en mathematique qu'en physique. 

Cote mathematique, il reste bien evidement a demontrer les asymptotiques, mais il reste aussi 
de nombreuses questions en suspens, comme la fonction tau isomonodromique. 

Cote physique, les matrices aleatoires ne sont qu'un outil. Le modele a deux matrices semble 
riche de perspectives pour mieux comprendre, voire pour donner un sens rigoureux, aux theories 
de champs conformes, en particulier sur les surfaces avec bords [S3]. Les operateurs de bords dans 
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les theories conformes sont encore assez mal connus, et le modele a deux matrices semble fournir 
un tres bon outil pour avancer sur ce sujet, surtout depuis que Ton sait calculer des observables 
contenant des traces mixtes. 

D'une fagon plus lointaine et ambitieuse, il semble que de nombreux concepts physiques et 
mathematiques sont en train de confluer: il s'agit de comprendre en profondeur la relation 
entre integrabilite(s) (au sens isomonodromique, au sens Yang-Baxter, au sens Ansatz de Be- 
hte), geometrie algebrique, combinatoire,... Mais ceci est plus un espoir qu'une realite a l'heure 
actuelle... 
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